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Many  standard  methods  of  test  for  concrete  and   concrete 
aggregates  contain  little  or  no  infonoation  that  will  give  the  user 
of  the  test  inethod  a  quantitative  Indication  of  how  vrell  individual 
test  results  of  a  series  ought  to  agree  with  each  other  or  hov  well  test 
results  from  different  laboratories,  on  supposedly  "identical" 
san^les,  ought  to  agree.  Such  a  statement  might  be  called  a  "precision 
statement «"  It  has  been  suggested  that  some  form  of  precision  state- 
ment  be  Included  in  ASTM  Test  Methods. 

In  a  recent  report  to  Subcommittee  Il-a  (Evaluation  of  Data) 
of  AS1M  CcHmaittee  C-9  (ik)*   it  was  suggested  that  AS^  designations 
be  divided  into  fotxr  classes  or  groups.  They  are  as  follows: 
Group  I   -  Methods  of  tests  for  xohlch  it  appears  possible 
to  obtain  a  measure  of  repeatability  and 
reproducibility  having  little  or  no  sample 
variance  either  by  repeating  the  test  on  the 
same  sample  or  by  making  synthetic  samples a 
Group  II  -  Methods  for  which  the  measure  of  repeatability  and 
reproducibility  will  necessarily  include  a 
component  of  variability  introduced  by  sampling. 
That  Is,  the  method  cannot  be  applied  a  number 
of  times  to  the  same  sample  and  it  does  not  appear 
that  synthetic  samples  can  be  made. 
Group  III  -  Methods  of  test  in  which  multiple  specimens  are 
required.  These  methods  might  possibly  be 


*  Numbers  In  parenthesis  refers  to  references  listed  in  the  bibliography. 


classified  with  those  of  Group  IX  but  some 
differences  should  be  recognized.. 

Group  IV  -  Specifications  and  miscellaneous  designations 
for  vhlch  no  statement  of  precision  is  needed. 

In  any  test  method  there  are  many  variables,  some  of  which 
are  inherent  In  the  method  itself  and  others  that  are  the  result 
of  outside  influences  Several  of  the  components  of  variance  that 
are  present  can  be  easily  noted.  There  may  be  within- laboratory 
variance  caused  by  differences  in  equipment,  operators,  etCo  There 
may  be  between-laboratory  variance  and  variance  caused  by  sample- 
to-sample  differences 0  Finally  there  is  the  variability  introduced 
by  the  method  itself  and  that  due  to  differences  In  materials » 

Not  every  ASTM  designation  requires  precision  statements  of 
the  type  under  considerationo  It  is  those  "<,<,,   methods-of-test 
specifications  wherein  apparatus  and  procedure  is  defined  for 
evaluating  a  magnitude  or  a  property  that  precision  statements  are 
neededo  Such  precision  statements  should  define,  within  stated 
confidence  limits,  the  maKlmum  and  the  normal  deviation  between 
test  results  that  may  be  expected  of  the  test  method  in  question 
when  performed  by  experienced  operators."  (^)>  The  ASTM  Manual 
on  Quality  Control  of  Material  (2)  contains  a  section  entitled 
"Presenting  -t  Limits  of  Uncertainty  of  an  Observed  Average."  The 
necessary  components  of  a  precision  statement  can  be  determined 
using  these  two  references  as  a  basis.  They  seem  to  be:  The  +  limits 
for  a  variable,  the  confidence  level  at  which  these  limits  apply, 
and  the  degrees  of  freedom  associated  with  the  determination  of  the 
limits.  Since  these  limits  are  determined  by  statistical  analysis 


of  Che  data  and  sample  statistics  are  used  It  is  accessary  to  use 
the  "Student's  t"  in  determining  confidence  limits o  This  is 
mentioned  by  Kaplan  (ll)  in  his  analysis  and  is  generally  discussed 
in  detail  in  most  statistical  texts  (ioOo  l^)o 

A  precision  statement  can  be  of  great  value  to  the  user  of 
the  test.  The  limits  provide  a  range  of  values  in  i^hich  the  results 
of  a  duplicate  test  should  lie  with  a  certain  probabilltyo  The 
operator  is  given  a  criterion  for  deciding  if  something  has  gone 
vrongo  The  variability  of  the  test  results,  vrhen  compared  with  the 
expected  variability,  can  give  an  indication  of  how  well  the  test 
is  being  performed.  That  is,  is  the  operator  using  the  proper 
techniques  and  is  he  being  accurate  enough  and  careful  enough  in 
his  measurements? 

If  the  repeatability  and  reproducibility  of  the  test  method 
is  kno^m,  it  is  then  possible  to  compare  results  of  different 
operators  and  laboratories  on  a  rational  basis.  Errors  that  are 
not  within  the  inherent  error  of  the  process  or  method  can  be 
discerned  and  with  a  kiiowledge  that  excessive  errors  are  present 
the  situation  can  be  analysed  accordingly. 


PURPOSE  AND  SCOPE 

The  primary  purpose  of  this  Investigation  was  to  Illustrate 
an  experimental  design  with  appropriate  analyses  t^  which  data  could 
be  collected  to  formulate  precision  statements  for  certain  classes  of 
Asm  Test  Methods.  A  secondary  purpose  or  a  result  deriving  from  the 
primary  purpose  was  the  development  of  precision  statements  for  the 
specific  methods  selected  for  Illustration.  Detailed  procedures 
are  given  which  Illustrate  the  general  approach. 

Statistical  analyses  of  the  data  were  made  and  from  these, 
confidence  limits  and  control  limits  were  determined  at  the  93   and  99 
percent  confidence  levels.  The  investigation  was  Intended  to  determine 
the  within- laboratory  precision  of  the  test  methods  for  different 
operators  following  exactly  the  procedures  for  the  methods  as  given 
in  the  ASIH  Standards.  During  the  testing  a  careful  analysis  of  the 
testing  procedure  was  made  to  determine  if  the  method  is  adequate  and 
if  the  desired  results  are  achieved  by  the  procedure. 

Since  precision  statements  (as  they  have  been  defined)  contain 
confidence  Intervals  and  degrees  of  freedom  it  was  necessary  to 
evaluate  the  means  and  variances  of  the  test  results.  From  these, 
linlts  on  the  sample  statistics  were  determined  for  the  two  confidence 
levels . 
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Three  ASIM  Test  Methods  falling  In  Group  I  (as  previously 
defined)  were  used,  namely; 


1.  C  117-^9  Test  for  Amount  of  Material  Finer  than 

No.  200  Sieve  In  Aggregate 

2.  C  127-59  7est  for  Specific  Gravity  and  Absorption 

of  Coarse  Aggregate 

3.  C  128-57  Test  for  Specific  Gravity  and  Absorption 

of  Fine  Aggregates. 


STATISTICAL  PROCEDURES 

To  accomplish  the  stated  objective  there  are,  generally  speaking, 
many  experimental  designs  that  could  be  used.  The  choice  of  the 
model  must  be  made  by  the  experlmentor  but  the  choice  will  necessarily 
be  influenced  by  the  Information  being  sought,  the  flexibility  and 
the  economy  of  the  model.  Careful  study  of  the  problem  before  testing 
has  begun  will  eliminate  the  gathering  of  data  that  will  be  useless 
for  the  testing  of  the  experimenter's  hypothesis. 

In  keeping  with  this  line  of  thought,  each  test  method  was 
examined  to  determine  variables  present  and  a  model  for  each  method 
was  selected  so  as  to  give  the  desired  information.  For  the  three 
methods  investigated,  the  use  of  factorial  models  seemed  to  be 
advantageous.  Ostle,  (12)  in  his  discussion  of  factorial  models 
states  the  following  advantages  and  disadvantages: 

Advantages: 

1.  Greater  efficiency  in  the  use  of  available  experimental 
resoiurces  is  achieved. 

2.  Information  is  obtained  about  the  various  interactions. 


3.  The  experimental  results  are  applicable  over  a  wider 
range  of  conditions;  that  is,  due  to  the  combining 
of  the  various  factors  in  one  experiment  the  results 
are  of  a  more  comprehensive  nature. 

k.     There  is  a  gain  due  to  the  hidden  replication  arising 
from  the  factorial  arrangement. 

Disadvantages : 

1.  The  experimental  setup  and  the  resulting  statistical 
analysis  are  more  complex. 

2.  With  a  large  number  of  treatment  combinations  the 
selection  of  homogeneous  experimental  units  becomes 
more  difficult. 

3.  Certain  of  the  treatment  combinations  may  be  of 
little  or  no  interest;  consequently,  some  of  the 
experimental  resources  may  be  wasted. 

It  appeared  that  the  advantages  far  outweigh  the  disadvantages* 
In  this  investigation  sampling  variation  was  eliminated  by 
performing  the  test  on  the  same  sample  a  number  of  times.  Since 
different  operators  would  be  using  the  test  method,  operators  were 
included  as  a  random  variable.  Type  of  material  was  also  Included 
as  a  random  variable.  In  addition,  the  amount  of  minus  200  material 
was  included  as  a  variable  for  test  method  C  117-'t^9*  The  range 
from  1^  to  lof)  minus  200  material  was  chosen  as  being  of  primary 
importance  and  four  levels  were  selected  within  this  range.  This 


variable  was  determined  to  be  a  fixed  variable  since  the  extremes 
were  represented  and  two  other  levels  covered  the  range  8u££lclentlyo 

The  following  is  the  mathematical  model  for  Test  Method  C  117-2|.9. 
The  mathematical  models  for  C  127-^9  and  C  128-57  are  similar  to  this 
model  except  that  there  are  t^;o  main  effects  rather  than  three o  This 
naturally  eliminates  all  terms  containing  the  main  effect  symbol. 
The  experimental  error  is  then  a  function  of  three  terms  (replicates, 
main  effect  1,  and  main  effect  2). 
Mathematical  Model  for  C  117-^9 

Y. ..  -  ■  an  individual  observation  in  the  ith  replicate,  in 

the  jth  level  of  factor  a,  in  the  kth  level  of  factor  b, 
and  in  the  1th  level  of  factor  c 
^   =  the  tjrue  mean  value 
P.  B  effect  of  the  ith  replicate 

0(<  ■■  effect  of  the  jth  level  of  factor  a  (material) 
R  *>■    effect  of  the  kth  level  of  factor  b  (amount) 
y.  ">  effect  of  the  1th  level  of  factor  c  (operators) 
(<^)..  a  effect  of  the  interaction  of  the  Jth  level  of  factor  a 

with  the  1th  level  of  factor  c 
((X^)..  "  effect  of  the  interaction  of  the  Jth  level  of  factor  a 

with  the  kth  level  of  factor  b 
^^^^kl  "  ^tt^ct.   of  the  interaction  of  the  kth  level  of  factor  b 
with  the  1th  level  of  factor  c 
(cx^^)...  >  effect  of  interaction  of  the  three  levels  of  factors 
^..  .  "  effect  of  the  azperlaiental  unit  in  the  1th  replicate 


8 
to  which  the  (jkl)th  treatment  combiixatlon  has  been 
randomly  assigned, 

and  the  terms  (Xit    }f -tf  i^o)ii   ^^  ^nkl  ^^^  assumed  to  be  independently 

2    2 
nonaally  distributed  with  escpected  values  of  0  and  variances  01  ,  (TL  , 

9 

2       2 

(^u,  and  0"  respectively. 


•^  i^^  •  t^-  'S^"^'^'=  -s^^^'"  -s<'"^'^j«  ■  ° 

Frior  to  testing,  an  Analysis  of  Variance  Table  (AN07)  was 
formulated  for  each  type  method  and  the  expected  mean  squares  (EMS) 
vrere  determined* 

A  Type  I  eanror,  usually  denoted  by  (K,  is  conmitted  if  the 
hypothesis  is  rejected  when  it  is  actually  true.  Obviously  one 
wishes  to  choose  a  test  in  which  the  probability  of  an  error  of 
this  kind  is  small.  It  has  been  found  generally  acceptable  in 
practice  to  use  an  (X -level  of  0.05>  althougli  this  depends  on  the 
type  of  problem  under  consideration.  For  this  st\uiy  an  (X  -level 
of  0.05  was  used  in  performing  F- tests. 

The  follo\7ing  equation  was  used  in  calculating  the  variance 
used  in  computing  the  control  limits. 


Variance  (Y,,,)  «  ^    +        ^+    *-+'"  + 


.   ,-      .       .    2  ^    2  ^   2 

Jkl    Ti^  n^     nj^p     n^     n^ 

where: 

2 
C^  B  Variance  component  due  to  o( 

2 

<T^   ■*  Variance  component  due  to  ^ 

2 
CC>  *"  Variance  component  due  to  f 

0^  "  Variance  component  due  to  the  interaction  of  oc  and  f 

2 
C^  c*  Variance  component  due  to  experimental  error 

and, 

n(^     "  nimiber  of  a's  ■•  a 

n  A      "  number  of  b*  s  ■•  b 


^aS      "  ^   times  b 

no   «  number  of  replicates  »  r 

n^   s  total  nuoobgr  of  observations  used  to  estimate 
variance  ('jfiij)*  ^  this  case  abr<> 

The  Variance  (^^1^1)  i-s  used  in  computing  control  limits  v/hlch 

are  placed  on  the  "grand'*  mean.  If  control  limits  are  placed  on 

means  which  are  averaged  over  only  on&   o£  the  variables  then  the 

appropriate  variance  (Y  )  must  be  used.  For  example;,  if  control 

limits  are  placed  on  means  averaged  over  replicates  only^  the 

(3-2    Q-2 

variance  would  be:  Variance  (Yp  )  •  — ^  +  — ~      .  l^ere  n 

\     n^     ng 

is  equal  to  the  number  of  observations  used  to  estimate  the  variance, 
in  this  case  n^  »n^«  ro 

PROCEDURES 

The  procedures  detailed  In  the  ASIM  Standards  were  followed 
exactly  by  the  operators  when  performing  the  tests  <>  Two  operators, 
acting  Independently,  performed  ten  replicates  for  each  test  method 
investigated.  The  operators  performed  several  practice  tests  to 
familiarize  themselves  with  the  method  and  test  procedures  before 
actual  testing  began. 

From  the  data  collected  by  the  operators,  sample  means  and 
sample  variances  were  calculated.  A  statistical  analysis  waa  made 
on  the  data  for  each  test  method.  Confidence  limits  were  determined 
for  9^  and  99  percent  confidence  levels.  From  the  statistical 
analysis  of  the  data,  precision  statements  of  various  kinds  could 
be  formulated  for  each  test  method  studied. 
Test  Method  C  117-^9 

The  purpose  of  C  117-lt'9  is  to  determine  the  amount  of  material 
finer  than  the  No.  200  sieve  slse  In  aggregate.  Briefly  the  test 
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consists  of  washing  a  tared  sample  of  aggregate  over  a  No.  200  sieve 
and  then  drying  It  to  a  constant  weight  and  determining,  by  weight, 
the  percentage  lost* 

The  material  used  In  the  performance  o£  this  test  was  local 
concrete  sand  of  glacial  origin,  de  sand  was  prepared  by  washing 
It  over  a  No.  200  sieve  several  times  to  elizalnate  all  inater:Lal  finer 
than  that  sieve  size.  The  sand  was  dried  to  a  constant  weight  and 
material  finer  than  the  No.  200  sieve  was  added  to  batches  of  it  in 
amounts  so  as  to  produce  ^00.0  gram  samples  containing  one,  four, 
seven,  and  ten  percent  material  finer  than  the  No.  200  sieve. 

Tt7o  types  of  minus  200  material  were  used  in  the  investigation. 
In  one  series  limestone  dust,  obtained  from  commercial  limestone 
filler,  was  used.  In  the  other  series,  a  limestone-residual  clay 
soil  was  the  source  of  the  fine  materials*  The  limestone  dust  was 
chosen  so  as  to  represent  a  non-cohesive  fine  material  and  the  clay 
was  chosen  so  as  to  represent  a  very  cohesive  fixte  material. 

For  each  replicate,  eight  samples  were  prepared  and  tested  in 
random  order.  Each  replicate  contained  samples  representing  the 
two  material  types  and  the  four  levels  of  each  material.  A  table 
of  random  numbers  was  used  to  achieve  the  random  order  <>  This 
precaution  v/as  taken  to  eliminate  any  bias  which  might  occur  due 
to  an  ordered  procedure  of  preparlog  and  testing  the  samples.  Replicates 
were  made  on  different  days  to  prevent  any  undue  Influence  from 
the  operator's  "mood"  on  any  one  particular  day  from  Influencing  the 
final  results. 


•  Ihiif led  Soil  Classification,  CH. 
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Test  Method  C  127-59 

The  purpose  of  C  127-59  is  to  determine  the  twenty-four  hour 
absorption  and  specific  gravity  of  coarse  aggregate «  Briefly  the 
test  consists  of  immersing  an  oven  dried  sample  of  coarse  aggregate 
In  water  for  twenty-four  hours,  then  determining  Its  weight  in  water 
and  its  weight  (SSD)  in  air.  The  sample  Is  then  dried  to  a  constant 
weight  and  ^^ighed.  The  absorption  and  specific  gravity  are  then 
calculated. 

TWO  materials  were  used  in  this  experiment o  One  was  a  crtished 
limestone  from  a  quarry  near  Bloomlngton,  Indiana,  the  other  was  a 
local  gravel o  The  coarse  aggregates  were  graded  so  as  to  have  equal 
amounts  of  material  retained  on  the  1,  3/^>  l/2ji  3/3  Inch,  and  Noo  ^i- 
sieves o  They  were  sieved,  washed  and  dried  to  a  constant  weight, 
then  the  various  fractions  were  combined  to  form  a  5000eO  gram  sasopleo 

For  each  replicate  the  two  sanq>les  were  tested  by  each  operators 
Replicate  determinations  were  made  on  different  days  to  prevent  biaSo 
Upon  completion  of  the  testing  the  t^^o  samples  were  again  sieved  and 
any  change  in  weight  of  each  fraction  was  noted. 

Test  Method  C  128-57 

The  purpose  of  C  128-57  is  t^o  determine  the  twenty-four  hour 
absorption  and  specific  gravity  of  fine  aggregate.  Briefly  the  method 
consists  of  soaking  a  tared  sample  for  tt^enty-four  hours  and  at  the 
end  of  this  period  drying  the  sample  to  a  saturated  surface  dry 
condition  where  upon  a  5OO0O  gram  sample  is  immediately  placed  in 
a  500  ml  flask 0  Enough  water  Is  added  to  bring  the  volume  to  5OO  ml 
and  the  weight  of  water  added  is  determined.  The  specific  gravities 
and  absorptions  are  calculated. 
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The  maCerial  used  for  chls  was  local  concrete  sand.  The  sand 
was  washed  to  remove  material  finer  than  the  Noo  200  sieve,  dried 
to  a  constant  weight,  and  separated  into  various  fractions  using  a 
logrithmlc  series  of  sieves <>  Sand  from  these  fractions  was  then 
recombined  to  form  samples  having  ti^ro  gradations,  one  with  a  fineness 
modulus  of  2oU0  and  the  other  a  fineness  modulus  of  3oOO<>  Each 
sample  weighed  1000.0  grams » 

Each  of  the  two  samples  was  then  tested  for  absorption  and 
specific  gravity  following  the  procedures  outlined  in  test  method 
C  128-^7«  One  procedtire  not  specifically  noted  in  the  ASTH  Standards 
was  used,  however o  l^en  the  sample  was  removed  from  the  ^00  mlo 
flask,  the  vrater  V7as  removed  by  decanting  the  sample «  The  sand  plus 
the  filter  paper  was  then  placed  In  the  drying  oven  together  and  the 
dry  weight  obtalnedo  The  filter  paper  was  oven  dried  prior  to  testing 
to  a  constant  weight  and  its  dry  weight  obtained^  By  subtracting 
this  weight  from  the  combined  dry  weight  the  sample  dry  weight  was 
obtained »  Tliis  procedure  was  used  to  speed  drying  time  and  while 
the  method  does  not  explicitly  state  that  the  procedure  should  be 
followed,  the  wording  does  not  preclude  its  use. 

The  portion  of  the  sand  not  used  in  the  flask  was  placed  in  a 
separate  container  and  dried  to  a  constant  weight  aloxig  with  the 
sample  from  the  flask.  After  the  dry  weights  were  determined  the 
portions  were  recombined  for  the  next  replicate »  As  before,  two 
operators  acting  independently  performed  ten  replicates  on  each  sample. 
One  hundred  percent  ethanol  was  used  in  minor  amounts  to  reduce  foaming 
in  the  flask. 
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Upon  completion  of  testing  the  txfo  saniples  v/ere  sieved  and  the 
change  In  weight  of  each  fraction  noted.  Loss  In  weight  occurred 
during  handling  of  the  sample > 

DISCUSSIOR  OF  RESULTS 

The  results  of  this  study  end  the  discussion  of  these  results 
have  been  divided  into  four  parts.  Tlie  first  part  Is  concerned  vlth 
the  analysis  of  the  data  for  Test  Method  C  117-49>  Tl^e  second  part 
is  concerned  with  the  axialysls  of  the  data  for  Test  Ifethod  C  127-59* 
the  third  part  with  the  analysis  of  data  for  Test  Method  C  120^31 > 
and  the  fourth  part  with  the  discussion  of  all  test  results  and  their 
application  to  the  development  of  precision  etatementso  The  second 
and  third  parts  have  been  further  divided  into  two  parts,  one  concerned 
with  the  analysis  of  absorption  data,  and  the  other  concerned  with 
the  analysis  of  specific  gravity  datao 

Results  of  Test  Method  C  117-^9 

It  should  be  recalled  that  a  S-^ACtor  factorial  model  was  chosen 
for  test  method  C  117-^9<>  Operators  and  material  were  random  variables 
and  amounts  of  minus  200  material  was  a  fixed  variable  <> 

A  ftartlett  test  for  homogeneity  of  variance  was  performed  on 
the  data  for  Operator  A  and  resulted  in  accepting  the  hypothesis 
that  the  within- sample  variances  were  equals  That  is,  the  eight 
variances  for  txra  materials  and  four  levels  of  each  were  accepted  as 
being  equal  at  an  C<-level  of  0.0^.  The  same  conclusion  was  reached 
for  the  within- sample  variances  of  Operator  B.  However,  a  Bartlett 
test  over  the  entire  table  of  data  (including  both  operators)  resulted 
in  the  rejection  of  the  hypothesis  that  all  wlthln-sample  variances 
rere  equal  at  the  0.0^  (X-level. 


Ik 


With  the  problem  of  aon-homogeneous  varlmice  present,  it  was 
deemed  advisable  to  analyze  the  results  of  Operator  A  and  Operator  B 
separately  since  one  of  the  assumptions  underlying  the  analysis  of 
variance  is  that  variances  are  homogeneous  ClO)o  Analyzing  the  data 
for  Operator  A  and  Operator  6  separately  eliminated  the  operator 
variable  and  hence  changed  the  statistical  model  from  a  3-£actor 
factorial  to  two  2-factor  factorial  models  similar  to  the  model  for 
analysis  of  test  method  C  127-59. 

The  calculations  for  the  sum  of  squares  for  the  analysis  of 
variance  tables  followed  the  procedure -presented  in  Ostle  (12)  page  351' 

Table  1  is  the  Analysis  of  Variance  Table  for  Operator  A,  Test 

o 
Hethod  C  117-^9»  The  estimate  of  experimental  error,  0"  ,  Is  ehawn 

to  be  0.01224^0  This  mean  square  is  used  in  testing  hypotheses  in 
regard  to  the  effect  of  replicates,  materials,  and  the  operator- 
material  interaction.  The  mean  square  for  amounts  was,  as  would  be 

expected,  very  large,  296.rrr3>  In  comparison  to  the  variance  t-Jhich 

2     .  % 

was  used  to  test  the  hypothesis  that  Go  »  0  (amount  effect }»  This 

effect  must  be  tested  using  the  mean  square  for  amount -material 

2 
interaction  since  the  expected  mean  square  for  amount  Includes  cr 

and  0^. 

F-tests  performed  on  the  data  from  the  analysis  of  variance 

indicated  that  materials  and  amounts  were  signiflcanto  That  Is  that 

2 
one  cannot,  at  the  0«05  oc-> level  accept  the  hypothesis  that  <7^  «  o 

(material  effect)  and  that  C7^  ■>  0  (amount  effect).  The  results  of 

the  F-tests  indicated  further  that  one  can  accept  the  hypothesis 

o  2 

that  (7^  »  0  (replicate  effect)  and  i^^a  ■"  0  (material -amount 

interaction)., 
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Table  2  presents  a.   similar  Analysis  of  Variance  Table  for 
Operator  B,  Test  Method  C  il7-49»  The  estimate  of  experimental 
error  was  0«0028o  Since  the  expected  mean  squares  for  this  design 
model  are  the  saise  as  those  for  the  previously  discussed  model  the 
Sisae  procedures  of  hypothesis  testing  were  followed.  As  before  the 

mean  square  for  amounts  ware  very  large,  298<'8232>  with  respect  to 

2 
the  mean  square  used  in  testing  the  hypothesis  that  C7^  »  0,  iiaxaely 

2     2 

1^"  teste  performed  on  the  data  from  the  analysis  of  variance 
of  the  data  for  Operator  B  indicated  that  amounts  xrere  significant, 
CJ^^  ^  0,  at  Oo05   (/-level o 

Amounts  would  be  expected  to  be  signi.ficant  since  their  means 
were  chosen  to  be  different  <>  The  fact  that  the  material  effect  was 
significant  for  Operator  A  icdleated  that  the  type  of  minus  200 
material  has  some  effect  on  the  results  obtained o  An  examination 
of  the  means  for  Operator  A  indicates  that  generally  a  higher  percentage 
of  minus  200  material  was  determined  for  the  finely  ground  limestone 
than  for  minus  200  clay  material  (Table  3)« 

TABI£  3 
UEANS  DETSmZNED  FOR  SAXA  OF  OPERATOR  A 


Material 

Limestone                                  Clay 

Percent  Added 

1             4          7           10           1           4          7           10 

Percent  Determined 
by  C  in-k9 

1.066  4,052  7-006    10.030  1.080  3,912  6,952    9-988 
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The  fact  that  material  was  found  not  significant  for  the  data 
of  Operator  B,  coupled  with  the  fact  that  the  variances  of  Operator  A 
and  Operator  B  were  found  to  be  different ^  indicates  that  the  results 
obtained  using  this  test  method  are  highly  sensitive  to  the  skill 
of  the  operator  performing  the  testo  Another  factor  lends  some  support 
to  this  suggestion.  Operator  B  was  an  "experienced"  operator  while 
Operator  A  X7as  not  (even  though  he  was  given  a  reasonable  "break- in" 
or  learning  period).  Therefore,  the  results  obtained  by  Operator  B 
may  be  indicative  of  results  obtained  when  the  test  method  is  performed 
by  an  experienced  laboratory  technician,  while  the  results  obtained 
by  Operator  A  may  be  indicative  of  results  obtained  when  the  test 
method  is  performed  by  a  person  not  trained  in  laboratory  techniques 
of  this  nat\ire» 

The  problem  of  establishing  an  acceptable  range  within  which 
results  of  duplicate  determination  should  agree,  or  establishing  an 
interval,  within  which  one  expects,  with  s«ae  given  probability, 
the  population  parameter  being  estimated  to  lie,  was  complicated  by 
the  non- homogeneity  of  variance  encountered* 

Another  factor  that  arises  is  select ii>g  the  type  of  limits  to 
be  set,  ioCo  confidence  llsilts  or  control  limits »  Confidence  limits 
on  a  mean  or  observation  provide  a  range  that,  with  a  given  probability, 
include  the  true  mean,  ^.  Control  limits  on  the  other  hand  provide 
limits  in  VThich  a  certain  percent  of  tloe  sanq>le  statistics  should 
fall  in  the  long  rtino  A  set  of  data  that  has  a  number  of  observations 
outside  the  control  limits  (greater  than  that  amount  determined  by 
choice  of  an  o(>level)  is  said  to  be  out  of  control.  Several  factors 
may  account  for  the  data  being  out  of  control.  One  can  be  the  fact 
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thac  some  asslsmablfi  cause  had  led  Co  discrepancies  bei:\7een 
observations,  but  another  is  that  the  precision  of  the  test  ixtethod 
or  process  is  incapable  of  producing  results  of  the  quality  required. 

For  the  purpose  of  this  study  It  seems  appropriate  to  present 
both  confldience  linits  and  control  limits »  lliese  limits  are 
presented  for  the  data  Operator  A  and  Operator  B  separatelyo 
Confidence  Limits  for  Operator  A 

s  «  OoOl2lt  from  Table  1,   therefore  s  «  Colli. 

Ninety- five  percent  confidence  limits  on  the  mean  of  n  observations 

made  to  determine  the  amount  of  material  finer  than  Noo  200  sieve  are 

s 
±  t^Q5  -=:,  wher^  t  is  based  on  the  chosen   o -level  (OoOJ  in  this 

case)  and  the  degrees  of  freedom  associated  x^ith  s«  In  this  case 

the  number  of  degrees  of  freedom  is  63 »  For  these  conditions,  the 

value  of  t  is  I.998  and  the  limits  are  Z"""^'-  « 

If,  for  example^  three  observations  are  made  the  confidence 
limits  placed  on  the  average  of  the  three  is  Y  t—^^-  »  Y  +  0«.128c 
The  interpretation  that  is  placed  on  these  limits  is  that  one  is 
9^  confident  that  the  true  mean,  |i,  lies  within  the  limits  established 
by,  Y  "t  O0I28  where  Y  is  the  observed  mean<> 

Control  Limits  for  Operator  A 

The  establishment  of  control  limits  for  a  measurement  process 
can  give  a  basis  for  determining  whether  observations  are  made  within 
the  limits  of  variability  allox/able  or  inherent  In  the  process.  The 
control  limits  were  placed  on  Individual  observationso  Since  both 
material  and  amount  were  determined  slgiiif  leant,  the  limits  were 
placed  about  the  average  of  ten  replicates «  The  limits,  placed  on 


20 


individual  observations  are  wider  than  limits  placed  on  means. 

The  limits  are  indicative  of  how  well  an  operator  luay  repeat  his 

measurenients  c 

The  limits  are  as  follows: 

Estimate  of  variance  =  <7v   °  -^    +  -211  »  0,013^ 

jr    1      1 

For  O0O5  a-level  based  on  n  observation,  T';!iere  n  equals  ten 

replicates,  ^  0«;  **  2*262.  The  control  li^aits  are  then  equal  to: 

t  2.262  sjO^Ol^    =  t  0.2623&. 

For  OoOl  (x-leveX,  t  «  3<'25  and  the  control  limits  are: 
+  3»25  \f~0o0l3k  «  t  0.376si» 

A  plot  of  the  observations  in  each  replicate  and  the  control 
limits  as  determined  above  are  shoim  in  Flgtore  1  for  the  clay  material. 
A  similar  chart  can  be  constructed  for  the  limestone  material  except 
that  the  mean,  X,  would  be  shifted » 

Confidence  Limits  for  Operator  B 

2 
s  »  0.0028  from  Table  2,  Therefore,  s  «  0.05^8. 

95^  confidence  limits  equal:  "t  ~fff=~  following  the  same  procedures 

shoxTT.  in  the  determination  of  confidence  limits  for  Operator  A. 
Control  Limits  for  Operator  "& 
Control  limits  for  data  of  Operator  B  vjere  placed  on  individual 
observations.  The  limits  were  plotted  about  means  averaged  over 
replicates  and  materials  since  amotmts  were  significant.  The 
procedures  follov/ed  in  determinixxg  these  limits  were  the  same  as 
those  used  for  determining  control  limits  for  the  data  of  Operator  A 
except  that,  obviously  a  different  variance  was  used,  that  associated 
with  Operator  B-   The  lisnltn   prer 
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*  0<ll6^  for  an    (x-level  of  O0O5 
t  0.159J&  for  an    <x-level  of  0«01 

Figure  2  shows  the  observations  and  the  control  limits  as 
determined  above* 

General  Observations 

Figure  1  and  Figure  2  give  an  Indication  of  how  wide  the  limits 
are  and  haw  some  of  the  data  were  scatteredo  This  scattering  could 
be  due  to  several  factors  (e.g.  a  variable  which  was  not  controlled^ 
such  as  laboratory  humidity,  may  be  signif leant) o  It  is  of  interest 
to  note  that  the  AS^  Coosslttee  responsible  for  this  test  method 
has  revised  It^  so  that  the  weight  determination  and  calculation  of 
percent  finer  than  No«  200  sieve  is  made  to  the  nearest  Ool  percent 
instead  of  the  nearest  0.02  percent »  The  remainder  of  the  revision 
appears  to  be  in  keeping  with  procedures  folloired  during  this  study, 
Zf  the  percent  finer  than  No.  2CX>  sieve  vras  determined  to  the  nearest 
0«1  percent,  the  variances  are,  for  practical  purposes,  aero  for  the 
data  of  Operator  B.  The  situation  is  then  similar  to  the  situation 
encountered  for  the  specific  gravity  data  of  Test  Method  C  l2T-59» 
Namely,  the  data  are  insensitive  to  statistical  analyses  and  prac- 
tical limits  of  1"  Ool  percent  seem  appropriate»  Noting  what  happened 
in  the  case  of  Operattnr  A  the  placeoient  of  these  "practical"  limits 
may  not  be  justified. 

Results  of  Test  Method  C  127-59 

Prior  to  testing,  the  variables  present  were  determined  and  a 
2-f«ctor  factorial  model  was  chosen  for  the  analysis  of  absorption 
and  specific  gravity  data. 

*  Accepted  as  a  tentative,  C  II7-6IT,  in  I96I. 
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Absorption 

A  BartletC  test  for  homogeneity  of  variance  was  perfosrnted  and 
resulted  in  accepting  the  hypothesis  that  within-sanple  variances 
were  equal  at  the  0«05  oC-levela  Ulth  these  variances  determined 
equal,  it  vaa   then  possible  to  proceed  with  the  analysis  of  the 
data  using  the  factorial  model. 

The  sums  of  squares  and  the  mean  squares  vere  calculated  according 
to  the  procedures  for  2- factor  factorial  models.  These  sisus  of 
squares  and  mean  squares  are  shoim  in  Table  4. 

F-tests  performed  on  the  data  from  the  analysis  of  variance 
indicated  that  there  x7as  &   significant  Interaction  between  materials 
and  operators,  i.e.  O^^a  ^  0.  Replicates  and  operators  were  not 

significant.  The  F~tests  also  indicated  that  materials  were  not 

2 

significant,  0^  «  0.  Materials  v/ere  espectsd  to  be  significant 

since  aggregates  having  different  absorptions  were  deliberately 
selected.  The  significant  interaction  can  be  logically  eKplalned. 
One  operator  may  be  more  experienced,  or  more  familiar,  with  one 
material  than  the  other  and  this  could  conceivably  influence  his 
results. 

With  estimates  of  the  variance  components  available,  limits  were 
determined^  As  before,  both  confidence  limits  and  control  limits 

were  determined*  Control  limits  were  placed  on  individual  observations. 

2 
The  calcul&tlon  of  the  components  of  variance  results  in  (T^  >  0. 

The  sum  of  squares  for  replicates  and  experimental  error  were  pooled 

to  form  a  better  estimate  of  the  population  variance.  This  variance 

was  determined  to  be  O.OOl'l-^  with  3^  degrees  of  freedom.  The  control 

limits  for  oc  0.05  and  (X*  0.01  are  indicated  below. 
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For  oc  «  0.0  5o 

0.05,  9  degrees  of  freedom  «»  2.262 

control  limits  «  t  2.262    J^^^^^^  "  -  O.O86S& 
For     OC"  O.Olo 

0*01^  9  degrees  of  freedom  »  3 '250 


control  limits  =  -  3.250     J    '  i     ^    „  +  0.133S& 


o 

Confidence  limits  were  determined  as  follows,  a  «  OoOOl^tJ 

(pooled),  therefore  a  »  O.O38.  For  this  case,  t  equals  2.029  based 
on  0.05  (X-level  and   36  degrees  of  freedomo  The  confidence  limits 
are  therefore,  t  2.029  ^^  =  *  %^  • 

Control  limits  were  also  determined  for  "between  operators." 
These  limits  were  determined  as  follows. 

The  estisoate  of  variance,  tT   ,  equals  0.000073*  ^e  limits  are 


-  t  \/0o 000073"  For  an  a-levol  of  0.0?  and  27  degrees  of  freedom, 
t  equals  2.052  and  the  limits  are  .  O.OiSsS.  For  sn    ot-level  of  0.01 
and  27  degrees  of  freedom,  t  equals  2.771  and  the  limits  are  +  0.02*^^. 

These  limits  are  narrower  than  the  control  limits  determined 
for  a  single  operator  and  are  applicable  to  situations  where  one  Is 
interested  in  making  comparisons  between  operator  means. 

Specific  Gravity 

A  statistical  analysis  of  the  specific  gravity  data  was  not 
warranted  since  the  range  in  values  for  each  type  of  specific  gravity 
never  exceeded  0.02  and  the  variances  were,  for  practical  purposes, 
zero.  The  practical  control  limits  for  this  case,  therefore,  are  +  0.01. 
These  limits  would  also  apply  to  between-operators  since  a  check  of 


27 

the  data  ladicated  almost  identical  specific  gravity  determinations 
by  each  operator, 

A  precision  statentent  as  it  mig^t  appear  in  the  test  toethod 
\70uld  be  as  follows:  Duplicate  detenainatlons  made  by  the  same, 
or  by  different  operators ,  should  lie  within  the  range  determined 
by  Y  I'OoOl,  where  Y  Is  the  observed  average  of  the  specific  gravity 
determinations. 

One   particular  point  was  noted  during  the  course  of  the  study 
of  Test  Method  C  127-59"  This  was  that  care  should  be  taken,  \ihea. 
submersing  the  sample  into  the  water  for  the  weight  determination 
in  water,  to  place  the  sample  in  the  water  In  such  a  manner  as  to 
prevent  as  much  as  possible  the  entrapment  of  aixo  Possibly  some 
reference  In  the  test  method  is  warranted* 

Results  of  Test  Method  C  128-57 

Prior  to  testing,  the  variables  present  were  determined  and 
a  2-factor  factorial  model  was  chosen  for  the  analysis  of  absorption 
and  specific  gravity  data.. 

Absorption 

A  Bartlett  test  for  homogeneity  of  variance  was  performed  on 
the  data«  The  test  resulted  in  accepting  the  hypothesis  that  within- 
sample  variances  xr&re   equal  at  the  0«05  o(- level  <>  IJith  these  variances 
determined  equal,  it  was  then  possible  to  proceed  with  the  analysis 
of  the  data  using  the  factorial  model o 

The  sums  of  squares  and  the  mean  squares  were  calculated  according 
to  the  procedures  for  2-factor  factorial  models.  These  sums  of  squares 
and  mean  squares  are  shown  in  Table  5<>  ^^  estimate  of  experimental 
error  is  0*022^8  as  determined  by  the  AMOV« 
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F- tests  performed  on  the  data  from  the  analysis  of  variance 

o 
Indicated  that  operator  effect  was  significant,  (Ja'  #  0,  and  also 

the  material  effect  was  significant,  C^,  #  0.  Replicate  effect  and 

the  Dperaftor -cater ial  interaction  were  determined  to  be  not  significant. 

Materials  might  be  eicpected  to  be  significant.  D3.fferent  fineness 
modulus  (different  gradations)  were  deliberately  chosen.  Since  the 
material  used  was  &  natural  sand  of  a  heterogeneous  nature,  it  is 
not  unusual  to  find  different  size  fractions  having  different  absorptions < 

It  is  interesting  to  note  that  operator  effect  was  significant, 
indicating  the  sensitivity  of  the  test  ciethod  to  this  variable «  The 
procedure  Involved  requires  judgsient  as  to  when  a  saturated  surface 
dry  condition  has  been  obtained  and  it  also  involves  sovm   amotint  of 
skill  on  the  part  of  the  operator  in  placing  500.0  granss  of  material 
into  the  flask  quickly.  It  should  be  noted  that  only  two  operators 
were  used  in  this  study  and  therefore  it  is  difficult  to  make 
generalities  about  the  whole  population  of  operators  with  confidence? 
^is  study  does  give  an  indication  as  to  what  might  be  expected. 

With  estimate  of  the  variances  components  ar/allable,  confidence 
and  control  limits  were  deterrained.  To  obtain  an  estimate  of  the 
precision  of  an  operator  the  control  limits  were  placed  on  the  individual 
observations*  The  procedures  followed  in  determining  the  control 
llnlts  vera  the  same  as  those  used  la  determinting  the  control  limits 

for  absorption  data  of  Test  Method  C  127-59*  For  this  case  a  pooled 

8  2  9 

estimate  of  a    was  determined  as  0.0190.  Since  C7^  -  0  and   O^p  «»  0 

these  were  used  with  a    to  determine  a  pooled  estimate  of  a  • 
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The  pooled  estimate  of  the  variance  Or  equals  0.0190  with  37 
degrees  of  freedom.  The  control  limits  were  determined  as  noted  belmr. 
For  ot=  0,05 

Control  liinits  »  I  2.^2  J^2^22  «  ±0.31^ 
For  »«=  0*01 

Control  limits  »  t  3.25  T^-^^^    «  t  C.kl^ 

Confidence  limits  were  determined  as  follcx^s.  The  pooled  estiiaates 
of  variance  equal  0.0190,  therefore  cr  equals  O.I38.  The  955^  confidence 
limits  are,  Y  *  2.028  —^5*^  »  Y  x  i^^  .  xhe  993&  confidence  limits 


are,  Yt  2.716  %^^  =  YI  ^j^I^. 
'  4n  >rn 


Control  iiaits  x-Tere  also  dsterminsd  for  "between  operators." 
These  limits  were  determined  calculating  the  variance  and  the  degrees 
of  freedom.  The  Satterwaithe  equation  was  used  to  determine  the 

degrees  of  freedom,  wliich  for  this  case  was  one.  The  control  limits 

+   r"2~       2 
were  determined  as  ^  t   la     •  The  a      for  this  esse  was  equal  to 

0.01325.  The  limits  were  then  1"  t  v' 0.01325  where  t  was  based  on  the 
chosen  oc -level  and  the  degrees  of  freedom  detsrained,  which  was  one. 
For  an  a-level  of  0.05  the  control  limits  are  t  I2.706  \f  6^013^  = 
.  I»k6/^*     These  limits  from  a  practical  viet^olnt  are  far  too  high. 
This  points  up  the  problems  associated  with  trying  to  make  generalities 
over  a  whole  population  with  only  one  degree  of  freedom  associated 
with  the  variance  component  due  to  operators.  It  should  be  further 
noted  that  the  operator  effect  was  determined  significant,  as  previously 
mentioned,  and  fact  would  tend  to  indicate  that  the  control  limits 
would  be  wider  for  this  case. 
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A  comparison  of  the  error  variances  of  this  test  isethod  with 
tliat  detenained  for  Test  Method  C  127-59  shows  that  this  test  method 
has  a  Hiuch  higlier  coiaponent  of  variance  due  to  experiioental  error 
than  did  Test  liathod  C  127-59»  These  error  variances  ^^ere  0*0226 
and  0.001^,  respectively. 

Specific  Gravity 

Results  of  i:he  armlyses  of  variance  perfonned  en  the  data  for 

the  three  types  of  specific  gravity  indicated  that  all  hiive  approsimately 

2 

the  saias  nagnitude  of  experinxintal  error,  nssely  a    =  0.0001.  HoT^ever, 

the  cstix'tSte  of  this  was  difficult  since  the  magnitude  of  computational 

figures  ere  siasll  and  the  factor  of  significant  figures  reduces  the 

2 
precision  of  the  analysis «  This  estiinate  of  a  however ,  if  used  in 

setting  control  limits  on  laeans  of  duplicate  detenaination,  replicates, 

if. 

gives  -  0,01*  This  Is  the  sasas  as  the  pr&ctical  limits  determined 

for  the  specific  gravities  of  coarse  aggregate  determined  by  Teat 
Method  C  127-59*  ^^e  specific  gravities  of  fine  aggrega?;©  as  determined 
by  Test  Method  C  128-57  ai^s  niore  variable  than  those  of  coarse  aggregates 
as  determined  by  Test  Method  C  127-59* 

Tlie  variances  associated  with  the  specific  gravity  data  were 
soiall  but  sosc  general  observations  may  be  made.  Bulk  specific  gravity 
had  more  variance  associated  with  It  and  is  more  sensitive  to  the 
operator  variable  and  to  the  material  variable  than  bulk  specific 
gravity  (SSO).  Bulk  gravity  (SSD)  in  tiim  is  sure  seiisltlve  to  the 
operator  variable  and  the  material  variable  than  is  apparent  specific 
gravity.  In  fact,  this  study  Indicated  no  variance  eOTq>onent  due  to 
the  operator  effect. 
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General  ObseinraCloas 

In  the  course  of  this  study,  it  was  noted  that  certain  factors 
in  Test  llethod  C  128-37  could  be  improved  or  elaborated  upon.  The 
teBq)erature  bath  specification  contains  no  provision  for  any  variance 
about  20  C.  The  time  specified  for  the  sample  to  recaain  in  the  bath 
also  has  no  tolerance  associated  x/lth  it.  The  flask  specified,  500  ml, 
is  actually  too  small  for  the  size  of  sample  being  testedo 

The  flask  Is  not  large  enough  to  allow  the  sand  to  be  mixed  by 
the  rolling  action  to  remove  entrapped  airo  It  vas  noted  that  it 
was  veiry  difficult  to  Tesaave   the  air  entrapped  in  the  send  at  the 
bottom  of  the  flask  due  to  the  fact  that  no  mixttire  vas  achieved 
by  the  rolling  action.  A  larger  flask  uould  provide  sufficient  room 
for  the  sample  to  be  mixed  x^ltli  the  water  and  the  entrapped  air  released » 

Ho  provision  Is  made  for  the  counter- act ion  of  foaming  that  may 
occur  during  the  removal  of  air  by  the  procedure  specified » 

Results  Applicable  to  the  Development 
of  Precision  Statements 
Both  confidence  limits  and  control  liioits  were  cosqputed  in  the 
statistical  analysis  of  the  test  data.  Control  limits  define  a  range  in 
values  vrlthln  which,  in  the  long  run,  the  sample  statistics  will  fall 
a  certain  percentage  of  the  tis^.  Confidence  limits  define  a  range  around 
the  sample  statistic  that,  with  a  given  probability,  will  include  the 
true  value  of  the  population  parameter.  Further,  ccmtrol  limits  ijidicate 

how  well  an  operator  is  performing  the  test*  If  the  operator's  results 
are  "out  of  control"  soma  assignable  cause  may  be  the  reason  (e.g.  the 
operator  may  not  be  following  exactly  the  test  procedures  specified). 


Si.  rtuo"  ©1' 


For  Inclusion  in  a  precision  stateoient^  control  limite  are  superior  to 
confidence  limits  since  they  provide  more  information  relative  to  the 
actual  performance  of  the  test  method. 

There  are  many  forms  which  a  precision  statement  may  take.  Ibe 
following  are  esec^les  of  several  types  that  seem  XTorthy  of  consideration 
for  use. 

1.  A  simple  statement  of  the  repeatability  and  reproducibility 

determined  for  the  test  method. 

Example:  Duplicate  determinations  should  check  to  the  nearest 
0.02  in  the  case  of  a  single  operator  and  to  the 
nearest  0.C5  for  between-operators. 

2.  A  statement  giving  the  estimate  of  experimental  error  detttrmined 

for  the  test  method^  the  estimate  of  each  variance  component, 

and  the  degrees  of  freedom  associated  with  each. 

Example:  The  esperictental  error  for  the  test  method  Is  0.0015 
with  27  degrees  of  freedom.  The  components  of  varl> 
ance  for  replicates^  operators,  and  materials  are 
0.001^  with  9  degrees  of  freedom,  O.OOS  with  1  degree 
of  freedom,  and  1.47  (highly  significant)  with  1 
degree  of  freedcm,  respectively. 

3*  A  statement  giving  the  estimate  of  esperimental  error  and  the 

eacpresslon  for  computing  control  limits  (e.g.  UCl  »  Y+tUQl  , 

IiCL  "  Y  -  t  JcSa    where  Y  =  the  sample  mean  and  Oo  "  the 

variance  about  the  mean).  The  degrees  of  freedc^  and  the  a- 

levels  to  be  used  wlili  t  must  also  be  noted. 

Example:  The  experimental  error  for  the  test  method  is  O.OOlg. 
Control  limits  are  to  be  computed  by  Y  ±  t  ^0~  , 
where  0^  "  <J£     on^  ^   is  based  on  an  « -level  of 
0.03  and  9  degrees  of  freedon. 

4.  A  statement  giving  the  control  limits  for  a  single  operator 

performing  the  test  a  number  of  times  and  the  control  limits 

for  between-operators .  Tine  degrees  of  freedom  and  the  (X -level 

for  determining  t  should  also  be  given. 
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Example:  Control  limlcs  for  &  single  operator  performing 

the  test  are  Y  t  0.03,  with  27  degrees  of  freedom 
associated  with  the  estimate  of  experimental 
error  axtd  t   based  on  9  degrees  of  freedom  and 
an  cC -level  of  0*0^.  Control  limits  for  betveen- 
operators  are  the  same  as  for  a  single  operator. 

Example  k  provides  the  optimum  information  and  seems  to  be 

the  best  for  inclusion  in  test  methods  of  the  type  considered  in 

this  study. 
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SUt&I&RY 


The  results  of  this  study  Indlcftte  that  it  is  essential  to 
detemlne  the  variables  for  a  test  method  and  to  develop  a 
statistical  model  before  actual  testing  begins »  The  estimate  of 
error  variance  and  ether  variance  components  can  be  determined  from 
an  analysis  of  variance »  With  this  infoinaation  available,  significance 
test  can  then  be  made«  Control  limits  can  be  placed  on  the  data 
of  a  single  operator  snd  on  data  for  bett7een->operator$o 

In  addition  it  was  found  that: 

a)  the  experimental  designs  used  produce  information  useful 
in  determining  control  limits  applicable  to  individual 
observations  of  the  attribute  being  measured* 

b)  control  limits  on  individual  observations  form  tbe  basis 
for  useful  precision  statements* 

c)  the  data  collected  are  also  useful  in  the  evaliuition  of  the 
test  method  itself. 

Precision  Statemento  for  Test  Methods  Studied 
The  following  are  cKamples  of  the  form  that  precision  statements 
mlgiht  take  using  the  results  of  thia  study*  These  statements  are 
based  on  the  previous  assumptions  and  include  components  of  variance 
due  to  replicates^  operators^  and  material*  ^fiiis  study  specifically 
eliminated  the  component  of  variance  due  to  sampling  and  other 
variables  not  mentioned*  The  precision  statements  are  for  test  data 
when  the  tests  are  performed  by  eaqwrieneed  personnel  following  exactly 
the  procedures  as  outlined  in  the  ASm  Standards  for  the  test  method. 
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Test  Method  C  117-2^9:  Test  for  >Iaterlal  Finer 
than  the  Ho.  SOO  Sieve  in  Aggregate 
Control  limits  for  a  single  operator  performing  the  test  are 
Y  t  0.12  percent,  with  63  degrees  of  tveedcm  associated  with  the 
estimate  of  experimental  error  and  t  based  on  9  degrees  of  freedtna 
and  an  oi -level  of  0*0^0 

Test  Method  C  127-59:  Test  for  Specific  Gravity 
and  Absorption  of  Coarse  Aggregate 
Absorption.  Control  limits  for  a  single  operator  performliig 
the  test  are  Y  t  0.09  percent,  with  27  degrees  of  freedcsa  associated 
with  the  estimate  of  experimental  error  and  t  based  on  9  degrees  of 
freedom  and  an  oC-level  of  0.05<>  Control  limits  for  between-operators 
are  the  same. 

Specific  C^avity.  Practical  control  limits  for  both  a  single 
operator  and  between-operators  are  Y  "t  0.01,  based  on  kO  observations. 
Test  Method  C  128-57:  Test  for  Specific  Gravity 
and  Absorption  of  Fine  Aggregate 
Absorption*  Control  limits  for  a  single  operator  performing 
the  test  are  Y  i  0.31  percent,  with  27  degrees  of  freedom  associated 
TTlth  the  estimate  of  e2q>erimental  error  and  t  based  on  9  degrees  of 
freedom  and  an  oC- level  of  0.05* 

Sx>ecific  Gravity.  Practical  control  limits  for  both  a  single 
operator  and  betmen-operators  are  Y  1"  0*01,  based  on  kO  observations. 
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Attached  is  an  informational  report  entitled  "Stress  and 
Deflection  in  an  Elastic  14ass  Under  Semiellipsoidal  Loads"  which  has 
been  prepared  by  VSto   John  L,  Sanborn  a  nieinber  of  the  staff  of  the 
School  of  Civil  Engineering.  The  research  reported  vas  the  thesis  of 
Mro  Sanborn  for  the  MSCE  degree  and  vas  performed  imder  the  direction 
of  Professor  E,  J,  Yoder. 

The  research  reported  is  concerned  ivith  the  estimation  of 
stresses  and  deflections  at  points  ifithin  a  pavement  due  to  applied 
wheel  loads.  Specifically  the  report  describes  the  application  of 
numerical  integration,  by  means  of  a  high  speed  digital  computer,  to 
the  solution  of  elastic  stress  and  deflection  equations  for  a 
semiellipsoidal  load  at  the  surface  of  a  uniform,  semi=-infinite  mass 
having  a  plane  boundary. 
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ABSTRACT 

Sanborn,  John  Leonard,  MSCE,  Purdue  University,  January  I965. 
Streae  and  Deflection  in  an  Elastic  Maes  Under  Semiellipsoidal  Loads. 
Major  Professor:  Eldon  J.  Yoder, 

A  common  problem  in  the  design  and  evaluation  of  highway  and 
airfield  pavements  is  the  estimation  of  theoretical  stresses  and 
deflections  at  points  within  the  pavement  due  to  applied  wheel  loads. 
Elastic  stress  distribution  in  a  homogeneous  mass  is  generally  used 
for  such  estimates,  but  its  application  in  this  field  has  been  limited 
to  assumptions  of  concentrated  load  or  of  uniformly  distributed  circular 
loads.  Formal  integration  of  the  expressions  for  stress  and  deflection 
for  non-unifoi^  loads  is  either  impossible  of  impracticable. 

Thie  report  describes  the  application  of  numerical  integration,  by 
means  of  a  high  speed  digital  computer,  to  the  solution  of  elastic 
stress  and  deflection  equations  for  a  semiellipsoidal  load  at  the 
surface  of  a  uniform,  semi-infinite  mass  having  a  plane  boundary.  Normal 
stresses  are  computed  by  the  Boussinesq  equations,  and  the  strain  due 
to  those  stresses  determined  by  accepted  elastic  theory.   Integration 
is  performed  by  Simpson's  Rule,  Three  sets  of  curves  representing 
stress  and  deflection  on  three  planes  normal  to  the  surface  are  included. 
The  curves  were  compiled  from  data  developed,  with  the  program  described, 
on  an  IBM  7090  computer. 

The  actual  program,  in  7090  Fortran  language,  and  an  example  of 
the  use  of  the  curves  are  appended  to  this  report. 
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INTRODUOTIOH 

One  of  the  most  frequently  encountered  probleme  In  design  and 
evaluation  of  flexible  highway  and  airfield  pavements  is  the  estimation 
of  stresses  and  deflections  at  points  within  the  structure  and  sub- 
grade  due  to  applied  wheel  loads.  Recent  studies  (19»  21)'*'  have  shown, 
in  particular, the  importance  of  determining  deflections  of  various 
components  of  pavements  in  order  to  evaluate  their  adequacy  or  service- 
ability.  In  the  past,  a  common  approach  to  this  problem  has  been  to 
consider  the  pavement  and  subgrade  to  act  as  an  elastic  mass  and  to 
evaluate  stresses  according  to  the  Boussinesq  theory  of  stress  dis- 
tribution. This  approach  implies,  of  course,  a  significant  simpli- 
fication of  the  layered  systems  actually  used  in  pavements.  The 
numerical  values  of  stress  and  deflection  determined  from  this  simpli- 
fication vary  significantly  from  true  values  both  as  a  result  of  the 
idealization  of  the  structure  and  of  the  difficulty  of  determining 
elastic  properties  of  the  materials.  Nevertheless,  patterns  of  measured 
stress  and  deflection  do  follow  theoretical  values  (21).  The  idealized 
case  may,  therefore,  be  used  reasonably  to  evaluate  relative  effects 
of  different  loads  or  different  pavements.   Basic  equations  for  stress 
due  to  a  concentrated  load  are  Integrated  over  the  loaded  area  to 
approximate  the  effect  of  a  wheel  load.  Strains  integrated  from  any 
given  depth  to  infinity  will  give  total  elastic  deflection  of  the 
given  point. 

The  major  difficulty  with  such  a  procedure  is  that  it  has  been 
*  Numbers  in  parentheses  refer  to  references  listed  in  Bibliography. 


practically  limited  to  uniform  circular  loads  because  of  the  unmanage- 
able expressions  vhich  it  is  otherwise  required  to  integrate.  A 
partial  solution  is  provided  by  influence  charts  (12,  15 )  prepared 
by  Newmark  based  on  the  effects  of  uniform  circular  loads.  These 
permit  determination  of  stresses  and  deflections  at  any  point  due  to 
any  uniform  load.  Even  this  requires  different  scale  drawings  for 
each  load  at  each  depth  for  which  stress  and/or  deflection  is  required, 
however.  Furthermore,  non-uniform  loads  require  superposition  of  a 
series  of  uniform  loads  to  approximate  the  effect  of  the  actual  load. 

Curves  of  stress  and  deflection  have  also  been  produced  (^), 
based  on  Newmark's  influence  charts,  which  are  practical  to  use  from 
a  design  or  investigation  standpoint.  These,  however,  serve  only  the 
uniform  circular  load. 

It  has  long  been  known  that  wheel  loads  vary  considerably  from 
circular  areas  of  uniform  pressure,  (7),  but  no  feasible  method  of 
evaluating  the  effect  of  these  loads  has  been  available.  With  the 
advent  of  high  speed  digital  computers,  however,  numerical  methods  of 
integrating  any  configuration  of  load  are  entirely  practicable.  An 
important  use  of  this  tool  is  to  compare  the  effects  of  known  load 
distributions  to  those  of  conventionally  assumed  distribxrtions,  knd 
to  provide  simple  and  useful  techniques  for  estimating  those  effects 
in  practice.  To  this  end,  the  computer  program  reported  here  was 
devised  to  evaluate  vertical  stress  and  deflection  according  to 
Boussinesq  theory  due  to  semiellipsoidal  loads.  The  completed 
program  was  used  to  prepare  curves  of  stress  and  deflection  which 
are  included  herein. 


PURPOSE  AND  score 

The  purpose  of  this  study  was  to  devise  a  numerical  solution  for 
theoretical  stresses  and  deflections  in  a  soil  mass  due  to  non-circular 
loads  of  non-uniform  distribution.  The  analysis  was  based  on  accepted 
theory  of  elasticity  and  previous  investigations  of  the  configuration 
of  ordinary  wheel  loads. 

The  study  was  composed  of  two  basic  parts— (1)  to  write  a  program 
for  a  digital  computer,  which  could  be  used  directly  or  with  minor 
modifications  for  any  given  load,  and  (2)  to  use  that  program  to 
develop  charts  of  vertical  stress  and  deflection  for  a  reasonable 
approximation  of  usual  tire  loads. 

The  computer  program  has  been  written  in  the  IBM  Fortran  language 
which  is  acceptable  to  a  variety  of  digital  computers  readily  available 
to  many  highway  and  airfield  agencies  as  well  as  to  most  universities 
and  colleges.  Variations  are  easily  introduced  to  provide  for  different 
load  conditions. 

Charts  of  stress  and  deflection  were  developed  for  a  semiellipsoidal 
load  of  proportions  which  approxi»at«  contact  pressures  and  areas  under 
some  standard  tires.  The  charts  are  in  the  form  of  curves  of  stress 
and  deflection,  respectively,  in  terms  of  maximum  contact  pressure, 
plotted  against  depth  in  teirme  of  size  of  contact  area.  Deflection  is 
treated  as  a  factor,  since  actual  deflections  depend  on  the  contact 
area  and  the  elastic  modulus  of  the  soil  as  well  as  the  contact  pressure. 

Stresses  are  considered  to  be  distributed  according  to  Boussinesq 


theory  which  assximea  a  homogeneoue.  Isotropic  elastic  mass*  No  attempt 
Is  made  to  account  for  the  effects  of  the  layered  structure  of  flexible 
pavements.  Several  investigators  have  described  solutions  to  two  and 
three  layer  problems.  However,  since  a  primary  function  of  this  program 
is  to  provide  a  basis  for  comparison  of  effects  of  non-uniform  loads 
with  the  conventionally  assumed  uniform  circular  loads,  the  yoA   vas 
restricted  to  the  simpler,  homogeneous  case  which  has  been  the  assumption 
adopted  by  other  inveetigators  in  a  large  portion  of  the  earlier  ifork* 
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REVIEW  OP  LITERATURE 

A  number  of  itnreeti£atlone  have  ehown  that  the  asetunption  of  elastic 
behavior  le  reasonable  for  the  response  of  flexible  pavements  to  traffic 
loads.  Foster  and  Fergus  (^)  compared  results  of  extensive  test 
measurements  on  a  clayey  silt  subgrade  to  theoretical  stresses  and 
deflections  and  reported  good  agreement*  Newmark's  Influence  charts 
(12,  1^)  have  been  used  extensively,  and  served  as  the  basis  for  Foster 
and  Ahlvln's  charts  (4)  of  stress  and  deflection  due  to  a  circular  load 
of  uniform  pressure  distribution.  The  latter  charts  are  costmonly 
used  In  design  and  research  problems  In  highway  and  airfield  pavements 
(20 }•  Palmer  and  Barber  (l4)  also  made  use  of  the  assumption  of  elastic 
behavior  with  satisfactory  results.  Baker  and  Papazlan  (1),  Walker  (19) 
and  Toder  (21 )  all  point  out  the  Important  contribution  of  the  elastic 
portion  of  pavement  deflections  In  determining  the  performance  of 
pavement. 

Burmister  (^)  has  considered  the  effect  of  the  layered  structure 
of  pavements  and  presented  curves  of  surface  deflection  at  the  center 
of  a  uniform  circular  load  for  the  two  layer  case.  Fox  (6)  extended 
the  work  of  Burmister  to  evaluate  stresses  at  various  points  In  the 
structure.  More  recently,  with  the  use  of  a  digital  computer,  Jones 
(9)  has  developed  tables  of  stresses  for  a  range  of  conditions  of  a 
three  layer  structure.  Tliese  are  all  based  oh  a  uniform  circular  load 
at  the  pavement  surface,  and  are  limited  to  points  on  the  axis  of 
Symmetry  of  the  load. 


There  have  been  relatively  few  atudlea  of  tire  contact  areae  and 
preaeurea.  However,  reaulta  of  several  atudlee  Indioate  elgnlf leant 
variations  from  the  aaauiaptlon  of  a  uniform  circular  load.  Among  the 
early  Investigator a  were  Heldt  (7)»  Teller  and  Buchanan  (17)  and 
McLeod  (11 )»  all  of  whom  note  the  nearly  elliptical  contact  area  of 
truck  tires.  More  recently,  Brahma  (2)  and  Lawton  (10)  have  made  the 
same  observations  with  respect  to  truck  and  aircraft  tires  respectively* 

Lawton' a  woxic  is  of  particular  intereat  in  that  both  contact  areas 
and  preaeurea  were  determizied  for  loads  on  a  aimulated  aubgrade  rather 
than  on  a  rigid  surface.  The  subgrade  used  was  a  mechanical  model  which 
has  been  deacribed  in  detail  by  Herner  and  Aldous  (8).  The  contact  areas 
for  a  variety  of  tlrea  and  loads  were  all  elliptical  and  had  a  reaaonably 
uniform  ratio  of  major  to  minor  axia,  averaging  about  1.7.  Preesure 
distribution  was  ellipsoidal  for  all  loads  within  rated  capacities  of 
the  tlrea,  and  varied  to  more  nearly  uniform  at  overloads. 

With  respect  to  numerical  methods  for  calculating  stresses  and 
deflections  in  an  elaatic  maae,  Stoll  (1^)  has  used  a  digital  computer 
for  calculation  of  vertical  stresses  in  a  foundation  soil  due  to 
variable  and  aaymDetric  column  loada.  His  results  indicated  that  the 
approach  is  aatiafactory  and  feaaible  for  use  in  similar  problems* 


PROCEDURE 

Computer  Program 
Program  Logic 
Th«  procedure  used  In  this  study  was  to  evaluate  streesee  in  an 
elastic  mass  by  numerical  integration  of  the  Bouseinesq  equations,  and 
to  compute  elastic  strains  resulting  from  those  stresses.  The  equations 
attributed  to  Boussinesq  express  stress  components  at  any  point  within 
a  semi-infinite,  elastic,  isotropic,  homogeneous  mass  due  to  a  single 
concentrated  load  on  the  plane  boundary  of  the  mass.  With  reference  to 
Pigxxre  1,  the  stress  components,  as  given  by  Taylor  (16),  are 

where  u   is  Poisson's  ratio  for  the  elastic  mass.  A  derivation  of  these 
equations  is  presented  by  Timoshenko  and  Ooodier  (16). 

The  elastic  strain  of  a  point  subjected  to  stress  is  determined 
from  the  normal  components  of  stress  by  the  equation 

where  z^  is  Poisson's  ratio  and  E  is  the  modulus  of  elasticity  of  the 
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FIGURE   I.      STRESSES  ON  AN  ELEMENT  DUE 
TO  A  CONCENTRATED  LOAD  AT 


THE  SURFACE 


material.  Total  deflection  of  a  given  point  dependa  not  only  on  the 
strain  of  that  point,  but  also  on  strains  of  other  points  deeper  in  the 
mass.  Therefore,  to  get  useful  deflections  it  is  necessary  to  integrate 
the  strains  on  a  vertical  line  from  the  given  point  to  infinity.   In  a 
numerical  procedure  the  suimnation  is  carried  out  to  such  depth  that 
the  remainder  is  negligible* 

In  applying  these  equations  to  distributed  loads,  the  basic  process 
is  to  express  the  stresses  and  strain  due  to  an  infinitesimal  element 
of  the  load,  treated  as  a  single  concentrated  load,  and  to  integrate 
those  expressions  over  the  entire  loaded  area.   It  is  evident  that  except 
for  symmetrical  loads  of  uniform  distribution  the  formal  integration  of 
the  expressions  is  extremely  andcward  and  impractical,  even  if  it  is 
possible,  for  a  given  load  configuration*  A  numerical  process,  however, 
using  a  high  speed  digital  computer,  is  entirely  feasible  for  any  load 
configuration,  the  only  requirement  being  to  define  the  area  of  load  and 
the  distribution  of  pressure  over  that  area.   In  this  process,  the  load 
is  divided  into  a  number  of  finite  elements,  each  element  being  treated 
as  a  concentrated  load,  and  the  numerical  integration  is  carried  out  in 
two  coordinate  directions  over  the  entire  loaded  area* 

The  numerical  integration  used  in  this  program  is  Simpson's  Rule  - 

S  =  t[yo"^  ^^^1"^  ^5  -^  •••  +^11-1^'*'  ^^^2"^  ^4  +  •••yn-2^-^  ^n]  ^^^ 
where  y^  are  values  of  the  function  at  equally  spaced  intervals,  h  is 
the  spacing  and  n  is  the  number  of  spaces  (must  be  even)* 

Figure  ^  illustrates  the  terminology  and  arrangement  used  in  develop- 
ment of  this  program.  Plane  XY  is  the  boundary  of  the  elastic  mass 
under  consideration*  Surface  ABODE  represents  the  distributed  load 
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applied  to  the  boundary,  and  Q,  a  point  in  the  maea  at  which  it  is 

required  to  determine  etress  and  deflection.  Coordinates  G  and  H 

locate  the  vertical  line  through  Q  with  respect  to  the  coordinate  axes 

XYZ,  and  point  Xjj,  Yjj  is  any  element  of  the  load,  Pjuv  represents 

the  maximum  ordinate  of  the  pressure  distribution  and  P^,  the  pressiure 

at  point  Xjj,  Yjj.  Rjj  is  the  horizontal  distance  from  the  vertical 

line  thro\i£h  Q  to  point  X..,  Yj„   and  is  defined  by  the  equation 

%  =  >/(Xij  -  H)2  +  (Yjj  -  0)2.  (7) 

The  program  as  described  here  provides  for  the  surface  ABCIS  to  be 

semiellipsoidal,  although  the  program  could  be  modified  to  permit  any 

desired  load  distribution.  Vertical  stress  is  determined  as  a  fraction 

of  P.,.„,  the  reference  preseure.  Deflection  is  treated  as  a  factor  P 
MAX 

In  the  equation 

A  '^~    P  (8) 

where  B  is  the  reference  dimension  and  E  is  the  modulus  of  elasticity 
of  the  soil.  Plow  of  the  entire  program  is  shown  in  block  diagram  form 
in  Pigure  J,  A  complete  Port ran  program  for  the  IBM  7090  is  given  in 
Appendix  A. 

The  first  step  in  the  program  is  to  map  the  pressure  distribution, 
P.J,  and  compute  the  radial  distances  R.j  for  each  element  of  load, 
and  store  these  quantities  in  two  arrays  in  computer  memory.  It  should 
be  noted  that  for  Simpson's  Rule  to  apply  an  even  number  of  equal  inter- 
vals, and  therefore  an  odd  number  of  points  in  the  loaded  area,  is 
necessary.  Computation  of  stresses  and  deflection  then  proceeds  along 
vertical  line  QS  (Pigure  2),  beginning  at  the  greatest  depth  used  and 
working  upward.  This  order  of  computation  is  used  in  order  to 
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FIGURE      2  A    DISTRIBUTED    LOAD 

ILLUSTRATING     TERMINOLOGY 
USED    IN    FORTRAN    PROGRAM 
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FIGURE    3.       FLOW     DIAGRAM 
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facilitate  integration  of  strains  and  to  eliminate  the  requirement  of 
storing  results  in  memory  until  a  vertical  line  is  completed* 

Because  the  orientation  of  the  point  relative  to  an  element  of  the 
load  is  variable,  it  is  necessary  to  compute  strain  directly  due  to 
each  element  of  the  load.  The  program,  therefore,  provides  for  computa- 
tion of  both  stresses  and  strain  at  a  single  point  of  interest  due  to 
a  single  element  of  load,  repeating  this  routine  for  each  element  of 
load,  then  sums  the  results  for  the  given  point  of  interest*  Similar 
computations  for  successive  points  along  line  QS  permit  the  determination 
of  total  deflections. 

Calculation  of  stresses  for  each  element  of  load  is  according  to 
equations  1,  2,  and  5*  Shear  stress  was  not  Included  since  it  has  no 
effect  on  vertical  strain,  though  it  could  be  included  easily  provided 
a  computer  with  adequate  memory  is  available.  With  normal  stresses 
known  it  is  possible  to  calculate  the  factor  <r,   -  xx^Tl -yu  (i^  in  equa- 
tion 5«  This  permits  the  factor  1/E  to  be  applied  Independently  for 
the  soil  conditions  pertaining  in  any  instance.  Because  the  strain 
required  is  that  in  a  vertical  direction,  the  factor  becomes 

Having  determined  stress  and  the  strain  factor  at  Q  due  to  one  element 
of  loadf  these  values  are  stored  in  computer  memory  and  the  routine  is 
repeated  for  the  next  element  in  the  array.  When  all  elements  in  one 
line  of  one  coordinate  dlrectlor  have  been  covered,  the  vertical  stress 
axid  the  strain  factor  are  Integrated  by  equation  6  and  the  result 
stored  in  memory.  The  same  process  is  used  tot   successive  lines  across 
the  loaded  area.  When  all  lines  of  load  elements  have  been  covered,  the 
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stress  and  strain  factor  of  the  total  load  is  calculated  by  another 
application  of  equation  6  in  the  second  coordinate  direction.  The 
total  stress  and  strain  factor  at  the  given  point  are  then  in  computer 
storage  available  for  printing  out.* 

Total  deflection  is,  as  noted,  the  stun  of  all  strains  along  the 
vertical  line  through  the  point  in  question.   Because  Simpson's  Rule 
requires  an  even  number  of  inteirvals,  it  is  not  possible  to  perform 
the  integration  at  every  point.  The  computation  of  deflection  factor  is, 
therefore,  limited  to  alternate  points  on  the  vertical  line.  Assuming 
zero  deflection  below  the  point  of  beginning,  an  initial  deflection  is 
calculated  for  the  first  point*  Depth,  stress,  strain  and  deflection 
factor  are  printed  out  for  this  point.  Depth  is  then  reduced  by  a 
predetermined  increment  and  the  program  returns  to  the  point  of  com- 
puting normal  stresses*  Vhen  stress  and  strain  for  this  point  are 
ccmputed,  they  are  printed  out,  with  no  deflection  computation,  and 
depth  is  incrsonented  again.  On  completing  the  third  run,  deflection  is 
again  computed,  taking  into  account  strains  of  all  points  computed 
previously,  and  again,  depth,  stress,  strain  and  deflection  factor  are 
printed  out.  This  process  is  repeated  continuously,  computing  deflection 
factor  only  at  alternate  points,  until  depth  is  essentially  zero.  No 
actual  result  can  be  calculated  for  depth  equal  to  zero  at  points 
beneath  the  load,  because  stress  expressions  for  this  case  become 
undefizMd  (i.e.,  zero  divided  by  zero).  End  points  for  the  stress 
curves  are  known,  however,  from  the  load  distribution* 

*  Horizontal  components  of  stress  were  not  stored  for  summation  because 
of  lack  of  memory  in  the  JEM  1620,  on  which  the  original  programming  was 
done*  It  should  be  noted,  however,  that  these  stresses  could  be  re- 
tained and  printed  out  with  larger  machines  such  as  the  IBM  7090  by 
adding  subscripts  to  their  notation,  reducing  to  a  fixed  orientation, 
and  repeating  the  summation  process  for  these  items* 
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This  completes  one  cycle  of  program  operation.  Program  control 
then  causes  reading  of  a  new  data  record  which  defines  a  new  vertical 
line  and  may  provide  new  load  configuration,  grid  size  and/or  Poisson's 
ratio.  The  computer  is  automatically  reset  to  initial  conditions  and 
the  entire  program  is  repeated.  Operation  continues  until  no  further 
data  are  supplied* 

Provision  is  made  in  the  program  for  a  series  of  three  different 
depth  increments.  This  is  desirable  because  the  curves  have  relatively 
sharp  curvature  near  the  surface  but  are  nearly  straight  at  considerable 
depth.  Using  large  increments  for  the  deep  points,  a  minimum  of  com- 
puter runs,  and  therfore  computer  time,  is  required.  In  areas  of 
sharper  curvature,  the  smaller  intervals  necessary  for  proper  definition 
of  the  curves  are  used. 

Input 

Two  input  formats  are  required  by  this  program.  One  provides 
data  relative  to  the  load  and  the  vertical  line  along  which  stresses 
and  deflection  factors  are  determined.  The  second  provides  the  depth 
boundaries  which  define  the  area  of  computation  and  the  depth  incre- 
ments in  each  interval. 

Data  in  the  first  input  record  arej  A,  the  ratio  of  major  to 
minor  axis  of  the  contact  area  (this  may  be  unnecessary,  or  inadequate, 
for  other  distributions  of  load);  N  and  M,  the  grid  size;  H  and  G,  the 
coordinates  which  locate  the  vertical  line  along  which  computation 
proceeds;  yUf  the  value  of  Poisson's  ration  to  be  used  in  the  computations; 
and  IDENT,  an  Identification  number  to  control  flow  of  the  program. 
The  second  input  record  Includes  DZl,  022,  DZ^,  ZPl,  ZP2,  and  ZF^  which 
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ar«  the  increments  of  depth  and  limiting  depths,  respectively,  for 
each  of  three  interval 8«  Computations  begin  at  ZF^,  end  depth  is 
incremented  by  DZJ.  When  depth  ZF2  is  reached,  the  increment  le  re- 
duced to  DZ2,  and  at  depth  ZFl  the  increment  is  reduced  again  to  DZl. 
The  use  of  two  formate  permits  keeping  input  to  a  minimvia.  Having 
once  established  the  depth  ranges  and  increments,  further  data  can  be 
proceeeed  without  redefining  these  values.  On  the  Initial  data  card 
a  negative  digit  is  punched  in  the  control  word  IDEMT*  This  digit 
causes  reading  of  the  following  card  which  must  contain  the  depth  data 
in  the  second  format  described.  Operation  then  proceeds  normally* 
In  successive  data  records,  as  long  as  the  depth  limits  are  to  remain 
unchanged,  the  identification  word  is  left  blank.  Program  operation 
then  by-passes  the  second  input  block,  proceeding  directly  with  com- 
putations for  the  new  line  of  interest*  New  depth  controls  can  be  in- 
troduced for  any  set  of  data  simply  by  punching  a  negative  digit  in  the 
identification  word  of  that  data  card  and  following  it  with  the  appro- 
priate depth  control  data  in  the  specified  format* 

All  linear  aeasurements  are  in  terms  of  a  characteristic  dimension 
of  the  loaded  area*  This  provides  dimensionless  results  for  compilation 
of  general  figures*  For  specific  locations  or  deflections,  however, 
results  must  be  multiplied  by  the  appropriate  dimension* 

Outpttt 
ha   each  set  of  input  data  is  read,  that  set  of  data  is  printed 
out  for  identification  of  the  results  which  follow  it.  Each  successive 
computer  run  throvigh  the  program  results  in  one  line  of  printed  output 
which  represents  conditions  at  a  single  point  of  interest*  I&formation 
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printed  Includes  depth  of  point  ae  well  as  vertical  stress  and  vertical 
strain  at  that  point.  Alternate  lines  also  include  the  total  deflection 
factor  for  the  point*  Thus,  complete  operation  on  one^data  set  pro- 
vides in  printed  output  the  coordinates  of  a  vertical  line  relative  to 
the  loaded  area,  and  a  series  of  results  for  points  along  that  line. 

When  all  points  defined  on  a  vertical  line  have  been  completed, 
the  program  automatically  reads  a  new  data  card  representing  a  new  line 
of  interest.  The  process  is  repeated  continuously  as  long  as  data  are 
supplied  to  the  computer* 

The  vertical  stress  indicated  in  the  output  represents  the  actual 
normal  component  of  stress,  at  the  point  described,  as  a  fraction  of 
the  majciawM  contact  pressure,  ^wav*  Deflections  are  expressed  as  a 
factor,  F,  which  oust  be  applied  in  equation  6  to  give  actvial  total 
deflections* 

Curves 
The  second  phase  of  this  study,  following  development  of  the  com- 
puter program,  was  to  develop  a  set  of  curves  of  vertical  stress  and  of 
deflection  factor  that  would  be  useful  both  for  design  application  and  for 
further  study.  A  number  of  difficulties  in  this  aim  are  immediately 
evident, among  the  most  important  of  which  are  - 

1,  Every  distribution  of  load  will  produce  a  different  family 
of  curves; 

2.  Asyiametric  loads  will  produce  a  unique  family  of  cxurves  for 
each  vertical  plane  through  any  reference  point; 

^*  Curves  of  deflection  will  depend  on  the  value  of  Foisson's 
ratio  used;  and 
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A,  Some  characteristic  pressure  must  be  known  as  a  base  for 

determining  actual  stress  and  deflection  for  any  given  load* 

There  is  nearly  general  agreement  among  investigators  that  tbe 
contact  area  voider  wheels  on  flexible  pavements  is  approximately 
elliptical.  Lavton's  work  on  aircraft  tires  (10)  confirms  this  and  also 
shows  that,  for  loads  within  the  rated  capacity  of  a  tire,  the  distribution 
of  pressure  is  essentially  semiellipsoidal.  The  proportions  of  the 
contact  area  are  fairly  constant,  the  ratio  of  major  to  minor  axis 
ranging  generally  from  1*6  to  l»d«  Maximum  pressure  depends  both  on 
tire  pressure  and  total  load. 

T^ical  curves  of  maximum  pressure,  ^vay*  ^b*  total  load  are 
shown  in  Figure  A,  It  can  be  seen  from  these  curves  that  at  rated 
load,  on  most  subgrades,  Pj^j^  is  practically  constant  at  1^2  per  cent 
of  the  average  rigid  plate  contact  pressure.  In  general,  the  average 
pressure  may  be  taken  as  the  inflation  pressure.  Thus,  P_^  would  be 
1.^2  times  the  Inflation  pressixre.  For  more  precision,  of  course,  it 
may  be  necessary  at  times  to  refer  to  actual  curves  of  contact  pressure 
vs.  inflation  pressxire  and  Piur/Pave  ^®*  ^^^^»   ^^   presented  by  Lawton. 

On  the  basis  of  available  information  it  was  decided  to  construct 
a  set  of  curves  for  a  semielllpsoidal  load,  using  a  ratio  of  major  to 
minor  axis  of  contact  area  equal  to  1.7*  This  is  shown  in  Figure  ^ 
where  A  =  1.7  snd  B  =  1.  The  equation  of  the  surface  which  defines  the 
pressure  distribution  is 

^  +  $+g   =1  (9) 

Hence,  the  pressure  at  any  point,  (I, J),  on  the  contact  area  ist 
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(a)     CONTACT   AREA 


SECTION     A-A  SECTION    B-B 

(b)    PRESSURE   DISTRIBUTION 

FIGURE     5.        CONTACT     AREA    AND    PRESSURE 

DISTRIBUTION 
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Ho  attempt  was  made  In  thia  etudy  to  evaluate  Polaaon'a  ratio  for 
any  particular  condition.  While  the  computer  program  providea  for  any 
yalue  between  0  and  0«^,  the  development  of  the  curvea  waa  in  conformance 
with  the  general  practice  of  taking  i^   =  0,5»  The  curvea  for  vertical 
atreaa  are  not  dependent  on  w  ,  of  courae,  ao  are  valid  for  any  aoil. 
The  curvea  for  deflection  factor  are  dependent  on^  ae  is  seen  in 
equations  2,  ^,  and  ^,  and  are,  therefore,  limited  to  the  assumption 
of  yU   s  0*^  which  implies  no  volume  change. 

The  question  of  orientation  is  not  as  easily  resolved.  Because  the 
load  is  not  distributed  symmetrically,  the  stress,  and  therefore  deflec- 
tion, a  given  distance  from  the  center  in  one  direction  is  different 
from  that  in  another  direction.  For  single  or  dual  wheels,  one  set 
of  cvirves  would  be  adequate  since  all  points  of  interest  would  be  under 
the  center  of  the  load  or  on  a  plane  through  the  transverse  axis.  When 
considering  tsmdem  axles,  however,  two  more  planes  are  immediately 
necessary  ->  one  along  the  longitudinal  axis  of  one  wheel  and  one 
throtigh  the  centers  of  two  diagonal  wheels.  Theoretically,  this  would 
require  a  different  diagonal  plane  for  each  geometrically  different 
arrangement  of  tandem  axlea  in  use.   It  would  be  possible  to  develop 
a  complete  set  of  charts  for  planea  oriented  eay  every  five  or  ten 
degrees  from  the  transverse  axis.  However,  since  one  wheel  of  a 
tandem  gear  has  a  relatively  small  effect  on  stresses  under  the 
diagonally  opposite  wheel,  the  result  would  seem  to  be  unwieldly  and 
confueing,  without  adding  any  great  usefulneaa  to  the  charts.  There- 
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fore,  one  diagonal  plana  was  selected  which  would  give  a  reasonable 
approximation  of  the  diagonal  for  standard  dual-tandem  gears.   It  is 
expected  that  for  most  practical  situations  the  actual  stresses  and 
deflections  determined  from  this  family  of  curves  will  give  values 
correct  within  the  precision  of  the  curves.  When  more  precise  work  is 
needed,  of  course,  the  computer  program  may  be  used  directly  or  a  new 
set  of  curves  may  be  produced  for  the  required  orientation. 

The  four  major  problems  are  thus  resolved  by  assuming  a  semi- 
ellipsoidal  load,  evaluating  curves  in  planes  oriented  in  three  directions, 
assuming  Poisson's  ration  -   0»^,  and  taking  the  reference  pressure 
as  the  maximum  of  the  ellipsoidal  distribution,  which  is  related  to 
the  inflation  pressvire. 

All  curves  are  kept  dimensionless  by  taking  the  characteristic 
dimension,  6,  as  unity*  By  this  technique,  the  curves  are  made  valid 
for  any  load  of  the  same  proportions  used  here,  regardless  of  size  of 
contact  area  or  magnitude  of  load*  Depths  and  offsets  measured  in  terms 
of  one  half  the  minor  axis  of  the  loaded  area  correspond  to  depths  and 
offsets  in  the  curves.  Results  are  plotted  for  depths  and  offsets  up 
to  ten  vinits.  Deflections  were  calculated  by  integrating  strains  to 
a  total  depth  of  5OO  units. 

The  completed  curves  appear  in  Figures  6  through  11.  The  orienta- 
tion of  the  plane  represented  by  each  chart  is  shown  on  the  figure. 
Figures  6  and  7  being  for  the  transverse  axis  of  load.  Figures  8  and 
9  for  the  diagonal  and  Figures  10  and  11  for  the  longitudinal  axis. 
Vertical  stresses  are  shown  in  Figures  6,  8  and  10,  Stress  in  per  cent 
of  maximum  contact  pressure,  Pj^j^,  ia  plotted  to  a  logarithmic  scale 
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on  the  horizontal  axis,  and  depth  to  an  arithmetic  scale  on  the  vertical 
axis.  Each  curve  represents  stress  at  a  distance  R  from  the  center  of 
the  load  along  the  line  indicated  on  the  figure.  Numbers  on  the  curves 
indicate  the  value  of  R  in  terms  of  B,  one  half  the  minor  axis  of  the 
contact  area. 

Figures  7»  9  and  11  present  deflection  factors  in  a  similar  manner. 
Actual  total  deflection  for  a  given  point  is  determined  by  equation  8. 
To  determine  stress  or  deflection  at  some  point  vinder  a  multiple 
arrangement  of  wheels,  the  appropriate  factors  are  read  from  the  curves 
for  each  wheel,  and  added.  The  sum  gives  the  effect  of  the  load  group. 
An  example  determination  of  stress  and  deflection  under  one  wheel  of 
a  tandem  axle  is  shown  in  Appendix  B* 
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Stn^MARY  AND  C0N0LU3I0N 

SuHmary 

This  study  has  produced  a  digital  computer  program  to  perform  a 
numerical  integration  for  determination  of  vertical  8tre8ee«  and 
deflections  in  pavements  under  wheel  loads. 

Loads  were  treated  as  semiellipsoidal  distributions  of  pressure 
over  the  tire  contact  area*   Bouseinesq  theory  of  stress  distribution 
was  applied,  using  Simpson's  rule  to  integrate  the  effect  of  the  entire 
loaded  area*  Elastic  strains  were  calculated  from  the  stresses,  and 
integrated  along  vertical  lines  to  determine  total  deflection  factors* 

Using  the  program  reported  here,  curves  of  stress  and  deflection 
factor  vs.  depth  were  plotted  for  three  different  orientations  of  a 
vertical  plane  through  the  center  of  the  load  -  transverse  and  longi- 
tudinal axee  and  one  diagonal  direction.  The  ratio  of  major  to  minor 
axis  of  load  was  taken  as  1.7*  The  maximum  pressure,  at  the  center  of 
the  load,  was  used  as  a  reference  pressure.  Poisson's  ratio  was  taken 
as  0*^*  Curves  are  provided  for  depths  and  offset  distances  up  to  five 
times  the  width  of  the  load* 

Re  c  ommend  at  i  ona 
The  following  recommendations  are  made  as  a  result  of  the  study 
reported  hereini  *: 

1*  That  further  study  be  made  of  the  distribution  of  pressure 
over  the  contact  areas  of  tires  -especially  of  truck  tires* 
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2.  That  the  ccmiputer  program  be  used  for  further  comparisons  of 
different  load  distributiona  -  in  particular,  the  effect  of 
uniform  loads  over  elliptic  contact  areas  compared  to  semi- 
ellipsoidal  loads. 

^,  That  a  study  of  theoretical  pavement  deflections  be  made 
comparing  results  of  three  layer  theory  with  those  of  this 
report* 
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APPENDIX  A 


Fortran  Program  for  IBM  7990  Computer 

DIMENSION  P2PI(25,25),R(25,25),U(25),V(25),W(25),U2(25).DY(25) 

110  FORMAT  (F5.2,2l4,2F7.2,F5.2,I5) 

111  FORMAT  (F5.2,5F7.2) 

112  FORMAT   (IH  ,F15.2,2F15.5) 
115     FORMAT   (IH  ,F15.2,5F15.5) 

5     READ  INPUT  TAPE  7,  HO,  A,N,M,H,G,RATIO,  IDENT 

IF(IDEKT)2,4,4 
2     READ  INPUT  TAPE  7,  111,DZ1,DZ2,DZ5,ZF1,  ZF2,  ZFJ 
A     WRITE  OUTPUT  TAPE  6, 110, A,N,M,H,G,RATIO 

DZ  =DZ5 

PDELT  «  0,0 

IffiLTA  =  0.0 

SUMD-O.O 

K  =  l 

IZ  =  1 

E  rM+1 

FrN+1 

DXr2.0*A/P 

X=DX-A 

DO  11  I'1,N 

TCMAX  =  3QRTP(1.0-X*X/(A*A)) 

DY(I)=  2.0*IMAXA 

TcDY(l)-YMAX 

DO    10    Jrl,M 

P2PI(I,J)  r(SQRTP(l,0-X*X/(A*A)-Y*Y))/(2.0»5.l4l6) 
R(I,J)  =  3QRTF((X-H)*(X-H)+(Y-G)*(Y-G)) 

10  YrY+DY(l) 

11  X  «  X+DX 
Ml  =  M-1 
Nl'  N-1 
Z  =  ZP5 

1     DO  16   Irl,N 
DO  15  J=1,M 

Qr  SQRTF(Z*Z*R(I,J)*R(I,J)) 
Q5  =  Q*Q*Q 
05  .  Q5*Q*Q 

3  r  R(  I,  J  )*R(  I,  J  >»-Z*Z+Z*Q 
V(J)  =5.0*P2PI(I,J)*Z*Z*Z/Q5 

V2  =  P2PI(I,J)*(5.0*R(I,J)*R(I,J)*Z/Q5-(1.0-2.0*RATI0)/S) 
V5^P2Pl(l,j)*(2.0*RATIO-1.0)»(Z/Q5-1.0/3) 
W(J) r  V(j)-RATI0*(V2+75) 
15     CONTINUE 
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I2rl 
69     EVENS  =  0.0 

DO  5  J=2,M1,2 

5  EVENS*EVENS*V(J) 
ODDS=  0.0 

DO  6  Jrl,M,2 

6  0DDS?0DDS4V(J) 

SUM  r  by(  I  )/5.0*(4 ,0*0DDS+2.0*EVEIIS  ) 

GO  TO  (72,75),  12 
72     U(I)-SUM 

I2s  2 

DO  7^  J'1,M 
Ik     V(J)  =  W(J) 

GO  TO  69 
75     U2(l)cSUM 
16     CONTHiUE 

15  =  1 

79  EVENS  =  0.0 
DO  7  I-2,N1,2 

7  EVENS  r  EVENS +0(1) 
ODDS  s-0.0 

DO  8  Irl,N,2 

8  0DDS^0DDS+U(I) 

SUM  r  DX/5  .  0*  (4 .  0*ODDS  +  2 .  0*EVENS  ) 

GO  TO  (77,78),  15 

77  SIGZ  =  SUM 
15  =  2 

DO  80  Irl,N 

80  U(I)tU2(I) 
GO  TO  79 

78  DDELT  =  SUM 

GO  TO  (50,51  ),K 

50  SUMD  =  SUMD+DDELT 

DELTA  =  SUMD*DZ/5.0+PDELT 

SUMD  =  SUMD+DDELT 

K -2 

WRITE  OUTPUT  TAPE  6, 115,Z,SIGZ,DDELT,DELTA 

GO  TO  (90,91,92),  IZ 

51  SUMD  =  SUMD4A.0*DDELT 

WRITE  OUTPUT  TAPE  6, 112,  Z, SIGZ,  DDELT 

K  =  l 

GO  TO  4l 

90  IP(Z-ZP2)100, 100,41 

100  DZ  »  DZ2 
PDELT  s  DELTA 
SUMD  -  DDELT 

101  IZ  *IZ+1 
4l  Z  »  Z-DZ 

GO  TO  1 

91  IP(Z-ZP1)1C2, 102,41 

102  DZ  =  DZl 
PDELT  a  DELTA 
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SUMD  =  DDELT 
GO  TO  101 
92     IF(Z-DZ)5,5,^1 
END 
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APPENDIX  B 

Example  Determination  of  Stress  and  Deflection 
Ab  an  example  of  the  use  of  the  cunres  presented  in  this  report, 
the  vertical  stress  and  deflection  at  a  depth  of  5^  in.  beneath  one 
wheel  of  a  dual-tandem  aircraft  gear,  at  rated  load,  is  determined. 
The  gear  configiuration,  as  shown  in  Figure  12,  consists  of  dual  wheels 
spaced  at  51-'~V^  in*  center  to  center  and  tandem  axles  at  61-1/4  in, 
center  to  center.  Total  gear  load  is  178,000  lb,  inflation  pressure  is 
176  psi,  and  contact  area  is  267  aq  in.  The  modulus  of  elasticity  of 
the  subgrade  is  2700psi. 

Contact  area  =  1.7TrB  =  267  aq  in. 


B 


-ff. 


267 


7TT 


=  7.07  in. 


^=  -^  -5.1 
B    7.07 


^MAX  "  ^'5^  ^  ^^^  ■  ^5^  P^^ 


Re  =  y  (51.25)2  +  (61.25)2  =  68.7  in. 


Wheel 

R 

R/B 

z(°/'Pmax) 

F 

A 

B 
0 
D 

0.0 
61.25 
68.7 
51.25 

0.0 
8.7 
9.7 
4.4 

6.00  (Fig.  6) 
0.25  (Fig.  10) 
0.15  (Fig.  8 J 
1.70  (Fig.  6) 

0.55  (Pig.  7) 
O.lt  (Fig.  11) 
0.0*-  (Fig.  9) 
0.19  (Fig.  7) 

2 

8.08 

0.62 

(fa  =  fg~  x  252  =  18.7  pei 


^  =  ^^Moo^^  0.62  =  0.58  in. 
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FIGURE    12.     EXAMPLE  GEAR  CONFIGURATION. 
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SYNOPSIS 


This  paper  reports  on  a  study  of  the  applicability  of  aerial 
photography  to  the  collection  and  analysis  of  hishtiray  traffic  floi-r  data© 
Specifically,  tiine-=lapse  aerial  photographs  i-rere  investigated  relative  to 
their  accuracy  and  efficiency  in  securing  speed,  volume,  headway  and 
other  pertinent  infornation  on  vehicle  iiB3ver.ient<, 

The  aerial  photographic  method,  althou^  iiapractical  for  the 
conventional  fixed  location,  limited  scope  traffic  survey  (e<,go  a 
volxuns  count  or  "spot"  speed  study),  demonstrated  a  significant 
potential  in  the  task  of  recording  a  wide  variety  of  traffic 
phenonena  over  "space"  as  vrell  as  tiiiBo  Time'-lapse  photography, 
■with  its  ability  to  record  a  driver's  behavior  in  a  series  of 
ei^qposures,  "was  Judged  to  be  the  most  accurate  and  practical  vjay  to 
analyze  vehicle  acceleration,  nerging,  divei^ging,  vreaving  and  passing 
patterns  over  an  extended  length  of  roadwayo  These  are  in^xjxtant 
traffic  flow  elements  i^Mch  are  all  but  injxjssible  to  study  by 
conventional  ground  methodso 

Acoong  the  nost  lii^x>rtant  assets  of  the  aerial  photographic 
method  is  its  ability  to  pictorial  ly  record  the  environment  in  vdiich 
traffic  data  are  obtainedo  This  attribute  enhances  the  value  of  the 
data  by  affording  the  investigator  idth  possible  reasons  for  unusiial 
traffic  behavioro 


AERIAL  PHOTOGRAPHY  APFLmU  TO  TR/VFFIC  STUDIES 

UlTHODUCTiai^' 

VJithin  the  next  decade^  laoro  thsxi  thirty  billion  dollaj's  t'dll 
be  spent  at  the  Fedextilj,  State  and  local  Isvels  for  the  construction 
of  nev!  roadways  in  tho  Liiited  States »     Additional  billions  xdJQ.  be 
allocated  for  the  operation  and  laaintciianee  of  existing  facilitioso 
Tovjard  \ftiat  ond  is  tliis  huo©  expenditure  boing  nade?    I-toat  highv;^ 
officials  \VDiild  agree  tl^iat  it  is  to  promote  the  efficient  flow  of 
the  nation's  goods  and  populations..    But  \Aiat  assurance  do  we  liav® 
that  this  program  \dIiJ.  actually  foster  efficient  transpcrtation? 
A  franic  reply  vjould  have  t©  b«^   "Very  little^  indeed"  o 

Piermrkably  little  is  kno'wn  about  the  patterns  of  flot^-  of 
motor  irahicles  on  the  highway  o    Tsnffic  data  on  speedy  voluiaB  and 
head'v-ray  eharacteristies  have  traditionally  been  gathored  at  a  fixed 
location^  thereby  denying  the  investigator  a  conplete  "Tiaw"  of 
flow  patterrjs  over  an  extended  length  of  roadirayc    Furthonnore^ 
raost  studies  have  been  concerned  vdth  but  one  or  tv»  t3:^fic 
chara<Jterlstics  and^  thus,*  fail  to  sreveal  possible  inter-relationshipa 
betvreen  the  raany  components  of  vehicle  flavio    Nuiaorous  attonqjts 
have  been  laade  in  recent  years  to  develop  theories  relating  the 
various  elements  of  traffic  inoveiQcnt  andj,  by  r.ieans  of  matheraatical 
nodels  aid  coir^iuter  teolmiquesj  to  apply  these  factoi*s  to  route 
capacity  detcr;ninatio\i  and  to  design  and  control  oritoriao     These 
effortSs  hoiievers  have  been  hampered  by  a  lack  of  appropriate  field 


dataj  ioOo,  data  collected  over  "space"  as  opposed  to  tl»  corrrontional 
"spot "9  or  £ixBd  looatlorif  studies  <> 

FoTjT  genei^  data  gathering  techniques  are  presently  In  uaso 
%■  far  the  most  conaaon  are  the  various  laanual  and  autoinatic  laethodsa 
These,  hovreTera  are  applicable  principally  to  "spot"  studies  of 
speed,  ToluEB  and  headvrays  and,  thus,  are  of  United  value  in  the 
study  of  vehicle  flow*  Jfovlng  test  vehicle  procedures,^  v.Mle  affording 
spaeial  data^  aire  quite  limited  in  the  accuracy  and  variety  of 
infonaation  obtainable » 

Photography  has  often  been  suggested  as  a  valuable  tool  in 
the  collection  of  spaeial  traffic  data©  As  early  as  1927,  the  speed, 
volune  and  headway  coE^wner^s  of  vehicle  flow  wsre  being  studied  vriLth 
the  aid  of  aerial  photography  (8)0  By  1935^  Bruce  Greonshields  had 
devoloj^ed  a  terrestrial  tina-lapse  motion  picture  technique  -vdiich 
yielded  acciu<ate  speeds  and  headv;ays  (4)0  Gresnshields  noted  in  1947 
the  potential  uses  of  aerial  photog3?aphy  in  traffic  research  (6)0 
In  recent  yearsg  Chicago  Aerial  Survey  has  pioneered  the  application 
of  stereo  continuous  stzip  photography  to  the  s'budy  of  traffic  flow 
(11,  9)0 

PURPOSE  A1\!D  SCOPE 
It  "uas  the  purpose  of  this  researeh  to  study,  in  tens  of 
applicability,  accural  and  efficiency,  the  tiioe-lapse  vertical  aerial 
l^otogj^aphy  nethod  of  oolleoting  traffic  flow  datao  The  technique 
1*08  investigated  relative  to  its  ability  to  detect  and  record  volune, 
speed,  vehicle  spacing  and  lane  use  infoxmatione  Consideration  vras 
also  given  its  use  in  gathering  density,  acceleration,  vehicle  placenent, 
nt^niffligrn  passing  distance,  and  other  pertinent  flow  datao 


rq   e/W- 


u^Mtiij  tov^^X 


DESCRIPTION  OF  THE  METHOD 
The  vertical  viox^  of  all  points  on  a  horizontal  plane  surfaoo 
yields  an  object  image  of  unchanging  scale  provided  a  single  height 
of  observation  above  the  surface  is  inaintainedo  The  aerial 
photography  technique  ponnlta  &  recording,  vjith  certain  inherent 
distortions,  of  this  image  on  a  photographic  filnio  The  scale  of  the 
photographic  iraage  is  a  function  of  the  lens  focal  length  (f )  and 
tlie  above  ground  elevation  (H)  of  the  camera  (aircraft)  as  equated  by 


Scale  «  s  »  H(feet) 
f (inches) 


Through  a  partial  overlapping  in  the  fields  of  view  of  tiro 
vertical  photos  takea  at  slightly  different^  but  laioitm,  points  iii 
tirae  and  sjsace,  the  motion  of  ground  objects,  such  as  motor  vehicles,) 
can  be  rieasured  and  recorded  in  terms  of  distance  per  unit  tirae  o 
In  addition^  measurenients  of  headways  between  moving  objects  and  their 
positions  relative  to  fixed  points  may  be  obtained  directly  froa 
the  aerial  photographs o 

Procedure  ♦»  D^t^  Gathezntn/^ 

The  Indiana  State  Highway  Conndsslon,  Bui-eau  of  Photogramraetric 
and  Electronic  Processes,  secured  and  processed  the  aerial  photography 
used  in  this  researcho  Eiqaloylng  a  K-iyC  aerial  camera  adapted  for 
a  high  speed  recycling  of  tvro  seconds  or  bettor  and  mounted  in  a 

State  oa^jned  Piper  Apaclie  aircraft,  nine  Inch  by  nine  inch  photograjiiy 
was  obtained  at  five  level  and  tangent  highway  locations  in  Haranond 
and  Indianapolis,  Indiana©  E|y  virtue  of  a  12-inch  focal  length 
lena  and  an  ISOO  foot  flight  elevaticaij,  a  photo  scale  of  approxlrnately 
150  feet  per  inch  was  malntoinodo 


As  ostimated  by  the  follovAng  relationship,  a  72  percent 
overlap  VQ3  considered  sufficient  to  assure  the  appearance  of  nearly 
all  vehicles  in  at  least  t;«5  photographs  •wiian  a  time  interval  of 
three  seconds  betvreen  exposures  T"jae  usedo 

Q.  iil*^<^^    (100) 
Is 

"tfiiere  Q,  =  percent  overlap 

1  "  the  photograph's  diniensionj,  in  inchess 
parallel  to  the  flight  line 

8  ■»  scale  of  photograjAy  in  feet  per  inch 

S6=  nvudnnun  expected  vehicle  speed  in  feet  jwr  second 

t  "  tints  interval  betvieen  eoq^osurss  in  seconds 

Settings        1  "  9  inches,  b  =  150  feet  per  inch^ 

S^''  100  feet  per  second  and  t  =  3  secojads, 
a  Q  of  72  percent  resultso 

This  inberval-overlap  combination  was  achieved  by  maintaining  an 

aircraft  ground  speed  of  approximately  S5  ii?>ho 

From  foio"  to  eight  tangent  flight  runs  of  one  to  two  miles 

in  length  were  nade  at  each  site  d\iilng  evening  periods  of  high  traffic 

volume  in  the  Sucua&r  ax^d  Fall  of  1962 »  By  flying  in  a  direction  opposite 

to  the  major  flow^  the  amount  of  data  obtained  per  run  \«as  ma^dmizedo 

Four  of  the  sites  were  recorded  on  Aerographio  black  and  v*iite  films 

The  Tri«State  Expressway,  a  four-lano  freeway  in  iwrthwestsm 
Indiana  (Figure  l)o 

The  Indiana  East^AJost  Toll  Road,  a  four-lane  trommj  in 
noirthvrestem  Indiana  (Figrrre  2), 
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Shadeland  AveniWi  a  lacjor  arterial  in  suburban  Indianapolis 
(Figure  3}o 

liadison  Avenue^  a  six-lans  e^qpresffMay  in  dovaatovai 
Indianapolis  (Figvura  4)o 

A  fii^h  location^  ^mdianapolis «  Ai-lington  AvenuBa  was  studied  using 

Aero-Ektachrome  color  film© 

The  data,  tiriie  of  day^  dxiration  of  noi^  shutter  speed  and  lena 
aperture  xrare  recorded  for  each  flight  pass*  A  tvK>  inan  crew  of  a 
pilot  and  photographer  proved  adequate  for  obtaining  the  desii'od 
datau 

Cttie  set  of  double  weight  contact  prints  v&s  made  for  use  in 
data  reduction  and  analysis o  A  second  set  of  identical  prints  on 
lightvraight  stock  vas  used  for  constructing  a  photo  index  of  each 
studj  sitao  Figure  5  illustrates,  at  a  reduced  scale,  the  standard 
nine  jjnch  by  nine  inch  photographs©  The  color  photography  vkls 
developed  as  a  roll  of  nine  inch  by  nine  inch  positive  transparenciaso 

j^'rooedure  ^^  Data  Induction 

Speed  and  distance  headway  qxxantities  vrere  obtained  from  the 
aerial  p^.otographs  by  sin^jle  linear  meaaurenents  en?)loyliig  a  100  unit 
per  inch  engineers^  scale o  A  hand  or  stand  niounbed  magnifying  glass 
vra,s  also  ;ised  by  some  of  the  data  processors  to  facilitate  the  precise 
nieasui^snjsnts  reqiiiredo 

The  black  and  white  photogi^aphs  were  fii-st  separated  into  sets 
representing  individual  fli^iht  passes  over  a  specific  highway  looationo 
Eaoh  exposiire  having  been  numbered  at  the  tirae  of  processings  the 
segregation  of  prJjits  was  rapidly  and  accurately  accon^ilisliedo 
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Kaoh  photo  set  vras  reduoad  independentlyo  First,  tho  scale 
of  photography  was  dotexroined  from  groimd  control  afforded  by  the 
paveaenb  a^qjansion  jolntso  J^  significant  change  in  scale  over  the 
length  )f  flight  ims  noted c  liarring  inajor  variationa,  the  mean 
scale  value  for  all  photos  of  a  given  site  and  day  \m.3   used  in  the 
data  analysiso 

YievriLng  the  photo  pairs  icLth  a  lens  stereoscope  pennitted 
an  apprjaiinately  level  section  of  tangent  liighway  to  be  accurately 
defined  for  each  siteo  The  time  duration  of  each  flight  run  having 
been  asiablished  in  the  field,  the  average  tioa  interval  betreen 
successive  esqjosures  \Jas   readily  eoii5)utedo  This  quantity  was  then 
average  i  over  all  the  flight  runs  iiade  at  a  given  site  and  intervalometer 
settings"^ 

The  vehicles  of  interest  -  those  flovdng  freely  on  the  through 
travel  lanes  ~  t/ere  identified  on  the  photographs  and  numbered 
successively^  by  direction  of  fiDivr,  from  the  start  of  the  photo  maxe 
Nuciberiig  \Ja5  done  directly  on  the  prints  vdlth  a  lead  or  grease 
penoilo 

Bach  "vehicle  appeared  in  at  least  tvjo  photographs o  Tliis 
perraitt 5d  its  speed  to  be  obtained  by  measuring  its  progress  along 
the  roadiv'ay  during  the  time  interval  between  photoso  A  location 
on  the  liighv^ay  coianon  to  both  photographs  served  as  a  reference  point 
to  i^Thioh  the  maasurenents  could  be  scaled  from  the  front  bun5)er  of  the 


*  The  ajstffi^lon  that  the  photographs  for  a  given  intervalometer  setting 
vssr©  squally  spaced  in  time  xfoa  verified  by  tests  conducted  on  tlie 
aiporift=<tanera=intervaloineter  system  under  simulated  fUglit  conditions< 
Variaiiojis  about  the  mean  interval  were  found  to  be  less  than  j*  0«>05 
seoondse 


vohiole  in  question  (see  Figure  6)0  The  pzocoduare  yielded  two 
values  which  "t«re  then  added  algebraically  30  as  to  give  the  distance 
traveled  (Dy)  in  the  photo  interval  o 

"HeadT-jay"  is  coniraonly  defined  in  traffic  engineering  as  being 
the  elapsed  time  bet^reen  the  passing  of  a  fixed  location  on  the 
highway  of  a  comtnon  point  (oogo  the  front  bun5)er)  on  tvx)  successive 
veliicleso  The  intenaitterrb  aerd.al  photographic  technique,  being 
essentially  a  moving  point  that  is  recording  vehicle  positions 
as  it  surveys  a  route^  cannot  yield  this  quantity  directiyc  Therefore, 
headvrays  were  initially  detennined  from  distance  neasureiiaents  as 
shovai  in  Figure  60  A  division  of  the  distance  headvray  by  the  speed 
of  the  "foUoxdng"  veliiolo  (Sf)  yielded  a  tirae  headway  (h-t)« 

■^■^ 

Traffic  density  (P)  tms  established  by  slnqjly  counting  the 
nuniber  of  vehicles  depicted  in  a  photograph  on  either  a  "by  lane" 
or  "directional"  basis o 

3inoe  each  axposurs  wasiaade  at  a  slightly  different  point 
in  time,  traffic  volumes  could  not  bo  detorrained  directly  from  a  count 
of  the  vehicles  picturedo  Hoi/ever,  a  single  relationship  between 
vehicle  and  aircraft  velocities  forms  the  basis  of  the  "Speed  Ratio" 
tochniqus  for  estimating  route  volumeso  This  procediire  considors 
volums  as  a  function  of  average  vehicle  speed  and  roxzte  densityo 

Utilizing  manual  or  automatic  ground  techniques,  volumes  are 
secured  by  recording  the  nunbor  of  vehicles  passing  a  fixed  location 
during  a  specified  period  of  tirae  o  Hovrever,  vAien  detennining  route 
voluEBs  from  aerial  photo ^raplis,  one  must  be  cognisant  of  the  fact 
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that  both  tha  riaferance  point  (or  aircraft)  and  the  vohioles  were 
in  ration*  Since  the  airborne  caiaera  did  not  photograph  tho  entire 
section  of  highvray,  A-B^  at  the  sane  Instant  in  time,  the  total  nuzaber  of 
distinct  vehicles  depicted  in  a  series  of  exposures  represented  not 
only  those  vehicles  in  the  study  section  at  the  flight  starts  biit 
also  those  v;hich  passed  point  B  over  the  dxiration  of  the  flight  run 
(see  Fifsure  d)o  It  -wls  this  voltuoe  at  B  during  a  tine  tp  vMoh 
the  Speed  Ratio  method  approximated  by  equation  1  i>1ien  the  aircraft 
•was  flying  opposite  the  direction  of  traffics 

^tp  "  Sp  ♦  S^      ('Hp)  o  o  o  «.  o    Bqo  1 

vdiere  V-^^  "  voluiae  of  traffic  passing  a 

^   fixed  point  in  time  tpo 

Sy>  "  average  vehicle  speed  in  miles  per  houro 

Sp  B  aircraft  ground  speed  in  adles  per  houro 

n^-  "  the  nuraber  of  vehicles  appearing  bet^veen 
^   principal  points  A  and  Bo 

Two  affSusqotions  vrere  inherent  in  the  development  of  equation  lo 
All  vehicles  "were  considered  to  be  traveling  ats  (l)  the  calculated 
average  speed  (Sy);  (2)  equal  distance  headiiayso 

If  the  aircraft  was  noving  directionally  vd.th  tho  traffic  but 
at  a  speed  greater  than  '3^^   the  voluns  of  traffic  passing  location  D 
in  tiua  tp  vaa  found  fron  the  relationships 


bauot 
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Time 


-^ >- 


^^=:Sd{        s, 


-^ >- 


C^-^C^  ^ 


Vehicle  per 
tp    seconds 


Study   Section    A-B 


^ ^-^ < ^ 


^^v 


Vt, 


•tp 


l^^  I^     , 


A-- Starting    point  of    flight    (photo   principal   point) 
B=end   point  of    flight    (photo    principal  point) 
BC-distonce    vehicles    travel   in  time    tp 
n.    =total  number   of   distinct  vehicles    appearing  on  the 
'P    olrphotos    between   A  ond   B    (vehicles  per  tp     seconds) 

Sp=  speed    of    plane  in   feet   per  second 

S^=  average   vehicle  speed   in   feet   per    second 

tp=  (t2-t,)--time  duration    for   flight    from   A  to   B  in  seconds 

t  =(t,-t,)  =  time,   in  seconds,   required    for   vehicle   n    °^    PO'"^  5 
'      ^     '      to  reach  point    B  when   traveling  at  S,.      Defined  as 

being    equal    to   tp. 
V»    --  vehicles    passing    point  B   in   time  tp 


FIGURE    7.     "SPEED    RATIO"  CONCEPT    FOR    TRAFFIC    VOLUMES 
FROM    STANDARD    AIRPHOTOS 

VEHICLES    AND    PLANE    MOVING    IN    OPPOSITE     DIRECTIONS 
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\  "     Sp  »  a^  ^"tp^     o   o   o  o   c      Eq.  2 

It  shoijld  be  noted  th-t  equation  2  is  aeaningloss  vdien  3p  <   Sy  o 
Data  reduction  from  the  Ektachrorae  aerial  iAotogra,phy  was 
acoon^jliahed  in  essentially  the  sane  manner  as  described  above  o 
In  this  oase,  hofr?errer,   positive  transparencies  were  vievred  on  a  light 
table  a  iUl  nuEiboring  and  maridjng    vas    done  vrlth  a  grease  pencil 
directly  on  the  "glossy"  side  of  the  photos » 

TRAFFIC  FlOV;  MUTIOIiSHIPS  STUDED 
i;iT«  Tr-S-L.\?Si;  .-iLliLU.  PlIOTOGIu\PPIY 

nearly  five  hundred  aerial  photographs  were  talcan  as  a  part 
of  the  present  researcho    The  nature  of  the  data  obtained  prohibited 
a  thorough  statistical  analysis  of  any  particular  traffic  characteristic 

or  hichvray  siteo    However,  wealth  of  information  renrealod  itself 

'a 

on  the  aerial  eaqwsttres,  indicating  the  potential  of  tixae^lapse  aerial 
photography  in  the  study  of  traffic  flowo  To  be  of  \xae  to  highway 
officials^  procedures  upost  be  capable  of  yielding  basic  information 
on  traffic  speeds^  headways  aiKi  volucieso 

Speede 

FiguTO  3  depicts  a  cumulative  frequency  graph  of  individual 
vehiola  speeds  conpllad  from  aerial  photographs  of  the  Tri-State 
fixpresmrayo  Th»  4ata  ww  gathered  over  a  one  mile  long  homogeneoufl 
section  of  the  highway  ao  as  to  afford  infoxmation  comparable  to  a 
"spot"  speed  study*  The  vehicle  speeds  of  Figure  8  were  plotted  for 
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several  diawotional  volunes  to  illustrate  the  tendency  for  velocities 
to  decrease  in  magnitude  and  range  i^ith  increasing  traffic » 

iiany  investigators  have  noted  that  average  speed  tends  to  vary 
as  a  function  of  the  traffic  volumeo  Ifean  speeds  for  various 
directional  volumes  of  traffic  on  the  Tri-State  Exparesavray  are 
shovjn  in  Figure  9s  along  vdth  a  best  fitting  straight  line  constructed 
by  the  method  of  least  squares  o  A  slight  decrease  in  speed  >jith 
increasing  volucsewas  indicatod© 

Time  Headwa?/Ta 

Cumulative  frequency  distributions  of  time  headways  for  tx*o 
lanes  of  uixLdlrectional  traffic  at  three  different  free  flowing  volumes 
■were  constructed  from  the  Tri=3tate  Espress^niy  photography  (see  Figiire  10) c 
A  comparison  of  these  curves  \vith  the  Polsson  distribution  seems  to 
confim  the  popular  belief  that  vehicles  are  randomly  distributed 
on  the  open  highirayo  As  expected,  the  frequency  of  short  headTvays 
increased  vdth  mounting  traffic  volvuaeso 

Lane  Distributions 

Vihen  studying  highvray  capacity  on  multilane  facilities^  the 
irrvostigator  is  often  interested  in  lane  usage  characteristics o 
Figure  11  relates  the  nunber  of  vehiolea  per  hour  traveling  in  the 
BBdian^  or  f&ssingy  lone  as  a  fraction  of  the  totals  tvfo-lane  directional 
volumso 
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The  p«ceding  disciussicn  has  been  confined  to  relationships 
involving  only  the  basic  traffic  flow  quantities  of  speed,  voIvum^ 
and  tin©  haadvjayo  Hovjever^  sevoxtil  investigators  (5^  7^  10^  et  n.lo) 
have  suggested  the  use  of  "spacial"  variables  as  being  more  meaningful 
and  accurate  representations  of  traffic  flow  phenononao  Data  for 
two  of  these  variables  —  traffic  density  and  acceleration—  are 
readily  obtained  with  tirae^lapse  aerial  photographyo 

Traffic  Density 

V/agner  and  IJay  (10)  have  siiggested  a  graphical  method  for 
pin-pointing,  in  terms  of  traffic  density,  locations  along  an  extended 
length  of  roadiiray  vihich  are  particularly  svisceptible  to  congestion^ 
and  to  determine,  in  units  of  time,  the  duration  of  these  hi^i  density 
conditionso  Utilizing  time^lapse  vertical  photography^ 
flight  runs  w©re  made  over  the  route  at  intervals  throughout  the  study 
periodo  Every  exposure  affoirded  a  density  for  each  direction  of  traffic  o 
These  densities liere  plotted  as  functions  of  time  and  distance  much  as 
one  \ould  construct  a  topographic  contour  mapo  Connecting  points 
of  equal  density,  Figure  12  depicts  contour  lines  at  intervals  of 
20  vehicles  per  lane  milso 

For  the  Ifedison  Avenue  example  illustrated,  liigh  densities, 
indicative  of  exbrem  congestiojn^  occxu7:9d  near  the  Pleasant  Run 
Paricway  signalized  interseotiono  This  condition  was  particularly 
severe  at  about  5:30  P.M.  when  vehicles  \»re  backed  vq?  nearly  1000 
feeto 
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FiSure  33  presents  average  speed  as  a  function  of  route  density 
for  eastbound  xaoveoents  on  the  Tri«Sto.te  Esqpressvjayo  Traffic  densities 
vrere  recorded  on  aerial  photographs  and  converted  into  vehicles  par 
direction  ndle  (all  lanes  in  one  direction)  for  analysiso  A  slight 
decrease  in  speed  \ib3  observed  to  aooompon^  itsoderate  incrooses  in  dorsityo 

Acceleration 

Althoijgh  traffio  density  is  a  function  of  both  speed  and  route 
volume,  it  alone  cannot  adequately  describe  the  nattire  of  vehicle  flowo 
At  aiQT  given  density^  the  facility  is  eithor  app2x>acliinci  con^ecrfcion 
or  energing  fi^>m  congestiono  This  dynamic  cor^xjsnent  of  traffic  flow- 
is  probably  best  described  by  vehicle  acceleration  pattcmso  A 
net  acceleration  of  vehicle  speeds  foretells  a  decrease  in  densi'ty, 
the  absence  of  acceleration  reveals  an  unchanging  density^ 
and  a  net  deeelex^tion  is  indicative  of  increasing  densityo 

The  tiiae-'lapse  aerial  photograjiiy  technique  is  an  ideal  means 
for  observing  the  acceleration  patterns  of  individual  motorists  as  they 
progress  along  the  highvray©  Figure  14  dlagracimatically  depicts  tli© 
variations  in  speed  for  each  of  five  vehicles  traveling  vrestbound 
on  the  Tri-State  Ejqiresswayo  Although  in  this  exan|)le  no  vehicle  was 
tracked  for  n»re  than  2500  f eet^  longer  traces  could  hove  been  achieved  by 
minimizing  the  speed  of  the  aircraft  relative  to  tlie  vehicle's  velocity 
ax^or  by  increasing  the  flight  elevaticaio 

Figure  15  relates  the  positions  of  the  five  vehicles  in  point 
of  ticB  as  irall  as  space o  Each  vehicle's  raovement  la  I'epresented 
by  a  line  vdth  the  slope  being  a  measure  of  spood  as  expressed  in  the 

equation. 

Tan  Q  "    ^distance        •  veliicle  s]^ed 

^  tilOB 
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FIGURE    15.     TRACKBJG  THE  lOVEl^iElIT  OF  ITJDIVIDUAL  VEHICLES  III  THE 
AIJD  SPACE -BY  THE  USE  OF  AERIAL  PHOTOGRAPHY 
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Changes  in  alope  (speed)  yield  veliiclo  accelarationg 


f^T^^      -l-^M^-^      «    vehicle  acceleration 
^  Time  >&  time 


lAnoa  cjTJSOinc  one  another  inr^doate  tliat  the  voliiclos  in  quoation  inare  at  the 
same  loiigitudinai  location  on  the  highvjay  at  a  concurrent  monent  in 
tiinsa    '>?hus#  either  tlra  vehicles  had   collided  or  onewita  passing  the 
otherc    The  plot  shoiirs  that  vehicle  nimiber  1  accelerated  from  ivl  mph 
to  56  raoh  in  passing  vehicle  n\ai4)er  3  oight  seconds  after  the  start 
of  the  ahoto  analysis  and  1100  feet  West  of  Korthcote  AvexiUBo    The 
passing  manouvsr  of  vehicle  nuniber  1  involved  its  moving  from  the 
outside  lane  to  the  median  lane  and  then  rstumingo    Tliis  vma  cooipleted 
in  appraxiiaately  13  seconds  and  in  sonBiihat  over  950  feerto 

Traffic  speeds  also  may  be  expressed  in  apacial  tonns  by 
plotting  the  average  speed  of  the  vehicles  appearing  in  each  eacposur© 
as  a  function  of  the  aerial  photograph  >s  location  on  the  route  o     This 
is  shoKn  in  Figure  16  for  the  Shadeland  Avenue  (State  Route  100)  site 
and  rsrveals  a  decided  increase  in  the  average  velocity  of  vehicles 
traveling  south  from  the  Pleasant  Run  Paxic.raj  intorsactiono    The  peak 
average  speadooearroddwit  10<X)  feet  b^nsnd  the  Parkway,  folloxired  by 
a  reduction  in  speed  as  the  drivers  approached  and  passed  l6th  3t*<eeto 
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Distance  Sr>acin/^ 

An  analysis  of  vehicle  spaciiig  patterns  assists  the  investigator 
in  deocribing  the  chaz^oter  of  traffic  density  andj,  as  such,  represents 
an  iii?x5rtanb  part  of  any  traffic  flovr  otudyo  Since  the  gap  (g) 
bet^-jeen  successive  vehicles  is  a  spaoial  quantity^  aerial  photography 
is  parti  culajrly  applicable  to  tliis  type  of  surveyo  Figure  17  depicts 
the  meaK  mdniraum  vehicle  spaoir^s  aUoxred  by  passenger  ear  raotoriste 
when  follo\iring  vehicles  at  various  speedso 

Figure  IB  5JLlustratos  the  x^sulta  of  a  study  of  lane  d-ianging 
as  a  function  of  directional  ti:^fic  volume  on  the  Tri=State  Expiressviayo 
The  data  vrere  obtained  from  aerieil  photographic  coverage  of  the  site 
by  noting  the  nuaber  of  vehicles  changing  lanes  as  a  fraction  of  the 
total  vehicle  count  for  each  flight  rario    A  regression  line  -vkuj 
statistically  fitted  to  a  few  scattered  points  and  indicated  an 
increase  in  the  frequency  of  lane  changing  voth  increasing  traffic 
flow  for  the  volumes  showio 


30 


175 


150 


/ 


-p 

<D 

O 


M 


•H 
O 

CO 
o 

•H 


.3 


125 


^ 


100 


75. 


S        50 


25 


□    1-fedian  Lane     observed 

^  .     . ,     T  Values ' 

o    Outside  Lane 


□ 


Passenger  Cars  Only 
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25  35  45  55  65  75 

Observed  Speed  (Sy)  of  Vehicle  Pairs  in  miles  per  hour 

FIGURE   17.     imJIMJl'I  SPACII^'GS  ALLOWED  BY  TH£  AVERAGE  DRIVER  WEN 
TRAlLIirc  AI-JOTHER  VEHICLE  AT  V/iRIGUS  SPEEDS  Oil  FOUR^ 
LANE  DIVIDED  HIGHWAYS  AS  f^ASURED  ON  AERIAL  PHOTO- 
GRAPHS 
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CSaJEEAL  iihUSlSlS  OP  TIE  ASRIAL  PHOTOGRAPHY  lEIHOD 
Data  Gathering  Characteristics 
Data  gathoriiTgj  as  the  tenn  is  engsloyed  here^,  enconmsses 
t^-jD  coii^wnents  of  ai^  traffic  surveying  method!  a  detecting  device 
and  the  recording  technique  o  It  is  toi^ard  an  analysis  of  the 
effioiencyy  accuracy  and  raeasurcment  precision  of  the  aerial  photo- 
graphic aothod  that  the  folloTdng  discussion  is  directedo  Although 
the  analysis  deals  prlioarily  with  the  colloction  of  speedy,  volume 
and  headiuay  data^  it  is  applicaile  to  ih©  study  of  all  traffic 
flow  characteristics  recorded  on  the  photograjihsa 

Equlpmsrit^,,  SugpJ^es  afxd  P^i^so.tyial 

The  previous  discussion  on  siethodolo^  has  indicated  the 
equipment,  SRipplies  and  pei-sonnel  angjloyod  in  collecting  data  by  the 
use  of  aerial  photographs  o  Those  items  are  summarized,  along  vrith 
estimated  costs,  in  Tables  1  and  2o 

Costs  may  be  esqiected  to  vary  somet^hat  idth  the  naturo  and  scope 
of  the  studyo  For  exarple,  the  costs  per  unit  of  data  gathered  v«)uld 
undoiibtedly  decrease  as  the  survey  period  is  lengthonedo 

Aerial  photography  can  be  enqsloyed  to  stuJy,  with  a  single 
con5)lement  of  equlpaient  and  personnel^  a  large  number  of  sites  within 
dose  proxiaiity  of  each  other  at  essentially  tho  some  tinso  Th« 
rasthod,  therefore^  becoues  increasingly  economical  as  the  study  locations 
beeomD  more  numerous  o  Color  photography,  however,  remains  significantly 
laore  e3q)ensiVB  than  the  standard  black  and  vAiite  aerial  photographyo 
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TABLE  1 

StiUIHiEIJT,  SUPPLDiS  ill©  PERSajlIEL  USED  II\I 

GATHERJl^G  TRAFFIC  DATA  «  BLACK  AND  VJHITB 

TUE-LAPSE  AZPJjlL  PHOTOGMPIII 


Basic  Eqtilpiaenbs 


Piper  Apache  PA«23  aircraft 

I&»17C  aerial  oaiaera 

12"  precision  lens 

B-3B  irrtervaloniBter  (2'«»120  seconds) 

ISarse  B»5  revriisd  type  film  processing  unit 


I&jor  Supplies? 


BzistJnan  Kodak  Super  XK  Aerographic  1623  black  and  vMte  film 

(or  Plux  X  iWrographic ) 
Film  processing  chemicala 
Photographic  paper 


Personn-els 


Aircraft  pilot 
Photographer 
Film  processor 
Photo  index  congjiler 


Unit  Costa  s 


Aircraft  (-with  caoera  system  and  pilot) 

Photographer 

Film  (150 «  roU) 

G175  e^qDOsures 

Film  processing  (labor,  chemicals,  etCo ) 
haavytreight  9"x9"  prints 
lightv/eight  9"3c9"  prints 


;y30oOOAour»* 
0  3oOOAour«fr 
i^57o75/rolL-»^ 
V  0o33/Q3q»3urB* 

U  0o26/print* 
U  0<.23/pnnt* 


Cotgjilation  of  index  sheet  (approxo  50  prints)  -      0  5«00/3heet"»«'- 


•*«•    Costs  estimated  by  the  Indiana  State  Highvray  Coraraission,  Bureau  of  Photo- 
grawnetric  and  Electronic  Processes,  lv62o 

i^i  Does  not  include  cost  of  the  9"x9"  prints  making  up  the  indexo 
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TABZiS  2 
EQUIPIMJT,   SUPPLIES  AiiD  PiJRSaillSL  USiiD 
IN  GATIIEimJG  TR/\FFIG  DATA  -  COLOR 
TliE-LAPSB  AERIAL  PHOTOGRAPHY  (TRAIISPAREIJCIES) 


Basic  I^quJLpissnbs 


Pipes?  Apacho  PA-23  aircraft 
IC=.17C  aerial  caiaera 
12"  precision  lens 
I^3B  intorvaloEJster  (2<=»120  seconds) 
Ibroe  B»5  remand  type  film  processing  unit  vith 
sevon  6  gallon  tanlcs  and  reel 


llajor  Supplies  $ 


Eastman  Kodak  Aero-Elcbaclironio  color  film 
CJolor  film  processing  chemicals  and  supplies 


Personnel! 


Aircraft  pilot 
Photographer 
Film  processor 


Unit  Costs? 


Aijroraft  (with  caaera  system  and  pilot) 

Photographer 

Fiaia(40«  roll) 

Q  45  exposures 
Film  processing  for  transparoncies 

Q  45  esqjosvures 


O30o00/hour«- 
si>  3o00/hour»* 
;i^72o00/roll^-- 
^  l<,6o/e3qjosuro 
!p46p90/roll-''^ 
iy  loQI»/o3q)osur8 


«    Costs  estimated  by  the  Indiana  State  Hlghv/ay  Commission,  Bureau  of 
Photograranetric  and  Electronic  Processes,  1962o 
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P^cj^lon  of  Recoxtled  U^tft 

For  -the  pxocedures  degoribed  in  earlier  sections  of  this 
report^  the  aerial  photographio  technique  afforcted  precisions  of  ^  Oo6 
feet  for  travel  distances,  ^o3  feet  for  distance  headways  and 
approximately  *0ol  seconds  for  the  photo  interval©  For  laseting  the 
reqtdra4iK»its  of  a  normal  traffic  study  these  levels  of  precision  are 
considei^  adequate  o 

Accuracy'-  of  Data  Gatherinx;  liethods 

The  aerial  photograjiiic  method^  >rfien  properly  exBcuted,  affords 
the  ♦i.O^  accuracy  required  of  nost  traffic  surveys  o  Knovjled^  of  the 
true  inlieiwal  betvreen  successive  photographs  and  a  precise  deteiraination 
of  the  scale  of  photography  are  fundacental  to  an  accurate  luoasurenianfc 
of  speedSj)  volumes  and  headways  from  tin»»lapse  |±iotogr^  hyo  In  the 
present  study,  the  interval  between  exposures  -was  approximately  three 
seconds^,  vdth  a  maximum  variance  over  a  given  flight  run  of  +00  05 
seconds  or  less  than  +2  percent  o 

The  scale  of  photography  "WBs  most  accurately  determined  by 
frequent  reference  to  established  groxind  control  points©  In  the 
present  studies  of  lefvel  higlrway  sectionss  the  enqjloysd  scale  ratio 
of  approximately  In  180(3^  once  precisely  d  efined,  vas  found  to  veiy 
less  tlian  *5  percent  over  the  one  to  two  mile  flight  run  or  across  aror 
individual  p^otographo 

The  nonaal  pitch  and  roll  effects  associated  with  standard 
aerial  photography  are  not  considered  significant  sources  of  err>©r 
in  scale  determinations  provided  the  pilot  is  ooq^orienced  in  aorial 
mapping  and  has  made  a  conscientious  effort  to  minimize  motion  about 
these  axaso  It  seems  probablttp  hovwverj  tliat  atmosirfioric  conditions 
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during  the  flight  would  liave  a  profound  effect  on  the  degree  of  scale 
variations  related  to  pitch,  roll  and  changes  in  flight  altitudeo 

Operational  Problems  and  IdiBi>ations 

During  the  course  of  the  research^  the  operational  problems 
and  limitations  associated  xdih  data  gathering  by  aerial  j^otogmphy 
viere  no-iiedo  Chief  aioong  tliese  viere  delays  due  to  inclenent  iieather^ 
the  method  8  s  restriction  to  daylight  hours^  obstructions  to  the  aerial 
view,  the  necessity  of  maintaining  a  tangent  flight  path,  and  a 
fluctuating  photo  scale  caused  by  changes  in  topograpl^  and/or  flight 
elevationo  In  addition,  va3?iations  in  plane  speed  occasionally  resulted 
in  insufficient  e^qposure  overlap^  the  critical  factor  in  assuring 
each  vehicle's  appearance  in  ttTO  successive  photographs o 

Data  Heduction  Characteristics 

The  reduction  and  tabtilation  of  data  is  usually  the  most  time 
consuming  aiKl  costly  phase  of  a  traffic  studyo  Thus,  it  is  appropriate 
that  soEB  attention  eftiould  be  devoted  to  an  analysis  of  the  economics, 
precision  and  tiras  efficiency  of  the  data  reduction  techniques 
invrestigatedo  For  discussion  purposes,  the  analysis  is  limited  to  the 
reduction  of  speed,  volume  and  headway  quantitioso 

EauiT^nt  and  Supplies 

Incbsigning  the  data  reduction  procedures  enqxLoyed  in  this  stuc^, 
every  effort  -^aa  oad*  to  minimize  the  amount  of  special  equipment  and 
skilled  labor  requiredo  Therefore,  the  methods  described  can  be 
readily  adopted  by  any  hlghxjay  organization  idLth  the  hiring  of  an 
adequate  clerical  staff o 


The  use  of  grease  pencils,  drafting  triangles,  appropriate 
engineers*  scales^  and  a  raagnifying  glass  brings  the  cost  of  reduction 
eqiiipraent  for  the  black  and  white  aerial  photographic  prints  to  about 
^•13000  per  mano    The  rsecessity  of  a  light  table  for  vievriLng  the 
transparencies  increases  the  coGt  to  ;^23oCX)  per  man  vrtien  obtaining 
data  frosi  color  exposuroso 

.Precision  of  Ifeasuraiasni^s 

The  precision  of  the  data  reduction  teolmique  associated  vdth 
the  photographic  raethod  •was  detciroiJiQd  in  light  of  the  requirements 
coaoon  in  traffic  analsrsis  and  the  capabilities  of  the  eqiiipcient 
engxLoyede    For  veliiclo  velocities  of  30  cdleo  por  hour  or  /^neator^ 
a  reduction  precision  of  at  least  ^  1  foot  per  second  or  j^  Oo7  milea 
per  hour    Ttfas      maintained  for  speed  data  obtained  from  150  foot  per 
inch  photographyo    LilcevJisei,  a  precision  of  ♦  Oo4  feet  was  realized 
for  dierbance  headways  and  other  linear  raaasurerientso    The  precision 
of  a  voluras  count  from  aerial  photograpliyvas  a  function  of  the 
time  duration  of  the  flight  as  vrell  as  systematic  errors  involved 
in  dorbeznining  vehicle  axKl  plane  speeds  o     It  is  estimated  that  the 
calculated  directional  voluraes  v»re  good  to  xdtlxin  ♦  10  per 
cent;,  vd.th  th«  precision  Improving  vrxth  increasing  flowso 

^•l^~^our  Becmirments 

A  coioplote  record  \ma  maintained  of  the  man-^hour  requLremente 
for  reducing  the  speedy,  volume  and  headvray  data  from  the  aoriel 
photog3ra.phSo     A  statistical  procedure  (tvro-way  classification,  mixed 
ncdel  analysis  of  variance)  vrais  erqployed  for  comi>aring  the  efficiencies 
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of  the  black  and  vAiita  VBirsus  the  color  aei-lalso     Tests  (F-test)  irer* 
raado^  at  a  O0O5  level  of  sicniTicanco,  of  the  hypothesia, 

Ua-Ub 

vdiere  JJj^  "  true  nwan  rate  of  data  reduction  from  black  and  vtfiite 
aerial  photographs*, 
Ufi  =  true  mean  rate  of  data  sroduction  from  the  color 
aerial  photogai^phso 

Vehicle  nxunberingj,  type  classification  and  velocity  measuremsnt 
rates  vare  combined  and  a  single  "speed  data"  reduction  rate  computed 
in  nindies  per  rehioleo  At  a  O0O5  level  of  probabilityj,  the  color 
and  the  black  and  >M.te  time^lapse  aerial  photography  nethods  vrere  not 
found  to  differ  significantly  in  their  efficiency  of  speed  data 
roductiono  The  mean  rates  for  the  black  and  white  photographs  and 
color  e3q»sxire3  vrere  2© 73  and  2o7h  minubes  per  vehicle  respectively© 

Testing  the  efficiency  of  the  headxray  neasuring  task  at  a 
O0O5  l€!Vel  of  significance^  the  black  and  idiite  and  the  color  photograplis 
differtfd  aignifioantly,  vdth  the  former  being  the  more  efficient  o 
The  roapeotive  reduction  rates  \mre  I0O5  and  lo30  minutes  per  vehicleo 

Variances  in  the  rates  of  volume  deteirainations  xvere  analyzed 
and  both  aerial  photographic  lasthods  vtera  found  to  have  an  average  rat« 
of  data  reduction  of  Oo53  minutes  per  vohicloo 

Pyolp3pr.m  and  Id^tationg 

liany  problems  unique  to  the  aorial  photographic  method  vrera 
observed  during  the  data  reduction  phase  of  the  research  0  In  addition 
to  "  Icaixig  "  vehicles  undor  bridges,  tho  pronoiuiced  tree  and  building 
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shadows  o£  the  late  affcomoon  often  ccJEplisated  VDhiolo  idontif ioation 
axii  iiapairtjd  the  pareciaioa  of  rasasxireanento  fjwm  the  photographso 
Similarly^  on  routes  running  East^lJeat  shadows  cast  by  the  vehiclea 
theiaselves  rendered  itieasuronisnts  to  the  daitened  bun^jera  difficult 
to  obtain,  Shadovwd  vehicles  vere  only  slightly  acre  distinct 
on  the  color  transparencies  than  on  black  and  white  printso 

Ibst  of  the  vrorloers  found  a  snail  magoifyinfj  glass  helpful 
in  identil^rinG  vehicles  and  measuring  distances  on  the  photos© 
Hovwver,  s«3(veral  coii5)lainsd  of  eye  stitiln  anl  attributed  it  to  optical 
distortions  in  the  lenso  Tiring  of  the  eyes  v/as  also  reported 
follovdng  a  prolonged  vievdug  of  the  color  transparencies  on  a  light 
table o 

AH  of  the  vrarlcers  admitted  being  distracted  by  the  roadside 
culture  depicted  on  the  photos  (eogo  a  large  outdoor  sidnming  pool 
adjacent  to  the  Tri-=Stat0  Expressway),  and  each  fo\aid  the  judgnente 
irequired  in  deciding  questionable  cases  of  vehicle  type  or  lane 
usage  to  be  time  consuiaingo  The  individual  black  and  irtiite  printe 
were  3omrvJha.t   ounibersono  to  work  vdLth  and  required  tedious  orientation 
before  neasurenents  could  be  obtainedo  Considerable  care  v*as  required 
vAien  measuring  distances  on  the  photographs  so  as  to  avoid  distortions 
due  to  parallax^  the  safest  pi-ocedure  being  to  laoacure  only  betvwen 
points  on  or  very  near  the  ground© 


Gaicuj'Jiou 

The  aeiial  photogrc?)hic  method  provBd  capable  of  affectively 
detoctinc  and  recording  the  basic  traffic  flow  elenonts  of  speedj, 
volmne  and  headvjayo 

In  obtaining  a  vdde  variety  of  additional  data^,  the  tin»- 
laps©  aerial  photographic  technique^  when  adopted  to  a  specific 
guanrey's  requireraentoj  vrould  appear  to  be  paarbiotilarly  lasefulo 
Such  physical  characteristics  as  vehicle  classification  and  roadvjay 
geoiaetjy  and  condition  are  readily  and  penaanenctly  recorded  l^  the 
photographic  methodo  Since  interrdttent  aeidal  photographs  afford 
a  ••viev;"  of  a  vehicle's  aovenent  over  both  space  and  tine,  the 
coEibined  physical  and  psyoholo-^ical  phenomena  in  traffic  flo\-r  may 
be  sttidied  in  terms  of  traffic  density}  acceleration  and  deceleration 
praotioesi  passing  behaviori  lateral  placementj  spacing  habitsj  and 
nergingp  divext^ing  and  reaving  patteznso 

AsaoiiQ  the  most  iEQArtant  assets  of  the  aerial  photographic 
nethod  is  its  abilSoy  to  pictorially  record  the  environment  in  vdiich 
traffic  data  is  obtainsdo  This  attribute  enhances  the  value  of  th« 
data  by  affordjjig  the  investigator  vdth  possible  reasons  for  unusual 
traffic  bobuvioro 

f^tandard  aericil  photography,  however,  caanot  be  efficiently 
enplo.vdd  for  a  "spot"  study  of  but  one  or  tvw  elements  of  traffic 
flo.r  at  a  fev/  sites  o  Photography  begins  to  have  practical  ^plication 
only  vAien  it  is  desired  to  procxares  concurrently,  a  permanent  record 
of  a  variety  of  traffic  and  roadt-ray  infoiroation  at  a  large  nu3.iber 
of  locations o 
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Tlia  use  of  standaxxi  aerial  photography  is  iirqiedeci  hj  the 
necessity  of  having  daylight^,  good  flying  ^■reathere  and  an  absence 
of  laajor  obatructlona  to  the  vertical  vieWo  In  addition,  the  coat 
of  collecting  and  reducing  basic  tjraiffic  data  on  aerial  photographs 
is  considerably  greater  than  vdth  conventional  giound  techniques 
In  all  but  the  moat  cojt5)lex  and  extensive  surveys  o 

Color  photography  yields  in  the  greatest  detail  a  coinplete 
vieiT  of  tlie  study  route,  its  vehicles  and  tlie  surrounding  culture  o 
The  advantages  of  color  over  black  and  tJhite  eoqjosures^,  however,  rarely 
Justify  tlie  additioxial  costo 

Aerial  photography  is  not  the  "ultimate "  in  techniques  for 
surveying  vehicle  traffic  characteristicso  In  cotmon  vdth  the  more 
conveiational  osthods,  photography  suffers  from  several  practical 
limitations  o  Wevortheloss,  the  techniqvie  affords  the  traffic  engineer 
a  valuable  tool  for  investigating  many  of  the  most  inportant  and 
coEQjlax  elements  of  traffic  flowo 
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